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Abstract 

This article starts by providing an introductory overview of the theoretical mechanics of 
rotating neutron stars as developped to account for the frequency variations, and particularly 
the discontinuous glitches, observed in pulsars. The theory suggests, and the observations 
seem to confirm, that an essential role is played by the interaction between the solid crust 
and inner layers whose superfluid nature allows them to rotate independently. However many 
significant details remain to be clarified, even in much studied cases such as the Crab and 
Vela. The second part of this article is more technical, concentrating on just one of the many 
physical aspects that needs further development, namely the provision of a satisfactorily 
relativistic (local but not microscopic) treatment of the effects of the neutron superfluidity 
that is involved. 



1 Elementary global mechanics of rotating neutron stars 



1.1 Introduction. 

Long before their observational detection as pulsars, theoreticiens were well aware 0] of the 
special physical interest of neutron stars - whose existence was confidently predicted - as 
well as of the (still entirely speculative) possibility of other more exotic (e.g. strange) stars 
of comparable compactness, meaning a radius only a few times larger than the Schwarzschild 
limit value, R = 2G M/c 2 , for a mass comparable with that of our Sun. Having presumably 
been formed by collapse of a stellar core that marginally exceeds the Chandrasekhar limit for 
for a self gravitating body with insufficient thermal pressure, a typical neutron star can be 
expected to a have a mass rather close to this limit, which - in terms of Newton's constant 
G, the speed of light c, the Dirac Planck constant H, and the proton mass m p - is given very 
roughly by the simple formula 

M~(!) 3/ V> (1) 

whose derivation is based just on the supposition that m p gives a rough estimate of the mass 
per cubic Fermi length, regardless whether the degenerate relativistic fermions in question are 
electrons (as in an ordinary white dwarf) neutrons, or even quarks. 

Unlike what was possible when superfluidity of the neutron matter in such compact stars 
was originally predicted [0] by Migdal, present day article accelerators can explore the physics 
of individual particle at energies that are now approaching the order of a TeV. Nevertheless, 
although their levels - from MeV to at most the order of GeV - are only moderate by such 
modern standards, the thermal energies - and particularly the Fermi energies - characteristic 
of matter in neutron stars remain beyond the range accessible in the laboratory for bulk 
matter. 

For a mass near the value given by (|l|), the condition that the stellar radius be large 
compared with the Schwarzschild value, R = 2GM/c 2 , places an upper bound 

P* « fl)V (2) 



on the mean stellar density - and hence also on the central density (since unlike what 
is possible other kinds of stars, a neutron star cannot have a density profile that is sharply 
peaked at the center). While less compact neutron star configurations (with lower mass and 
larger radius) can exist in principle, it is hard to see how they could be created in nature, 
so a typical example can be expected to have a central density that is not so very far below 
what is permitted by this Oppenheimer - Volkoff bound (|2|). Since this bound is interpretable 
as the order of a proton mass per cubic proton Compton length, it is evidently quite a lot 
higher than the density of the order of a proton mass per cubic pion Compton length that 
characterises ordinary nuclear matter. In terms of the pion mass m n this ordinary nuclear 
density will be given in order of magnitude by 

P nuc « (^) 3 m^m p , (3) 

which is a few times 10 gm/cm 3 . The prediction that typical neutron star core densities 
are thus well beyond what is easily accessible to experiment is one of the reasons why it is 
so interesting, not just for astronomy, but also for the basic physics M, H of bulk matter at 
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the corresponding intermediate energy levels, to acquire and analyse as much observational 
information as possible about neutron stars (as well as "strange" or other comparably compact 
stars, which, if they exist, will also have core densities in the same range). 

In addition to the limited amount of such information that is available from other mecha- 
nisms (such as binary orbital behaviour), we are fortunate to have at our disposal an enormous 
and steadily increasing body of relevant information provided ( see Figure |l]) by pulsar timing 
measurements: radio (and in some cases optical or other) observations provide continuous 
high (sometimes within 10 -9 ) precision monitoring of pulsar frequencies, which are generally 
believed to correspond directly to the rotation frequency f2 of the underlying star, or more 
precisely to that of its rigidly rotating outer "crust" . 

The present article (like a briefer preceeding review |5[) is intended as a self contained 
introduction to the theory of the phenomena most relevant to such observations. It is meant 
to be accessible to non-specialist readers, who are assumed just to have a grounding in gen- 
eral physics, at the level provided by Landau and Lifshitz in areas including relativistic 
gravitation theory, and non relativistic superfluidity and superconductivity theory. 

As discussed in detail in accompanying articles in this volume, outside a still mysteri- 
ous core (that may consist of quark matter) neutron stars are generally believed to consist 
essentially of a neutron fluid interior and a surrounding crust. The outer crust material is 
qualitatively similar to an ordinary metal, consisting of baryons concentrated (as a majority 
of neutrons with a minority of protons) in nuclei in a degenerate Fermi type sea of elec- 
trons at concentrations up to and beyond the white dwarf limit, where the electrons become 
relativistic, at a density given in terms of the electron mass m e by 

Pre\ « (^J rn^m p , (4) 

whose value, of the order of 10 7 gm/cm 3 , corresponds to about one proton mass per cubic 
electron Compton length. 

The transition to the qualitatively different kind of material that makes the behaviour 
of neutron stars so very different from that of ordinary degenerate electron supported white 
dwarf stars occurs at a critical "neutron drip" density pdvip that is reached when the Fermi 
energy of the degenerate electrons becomes comparable the binding energy -E nuc per baryon 
in a nucleus, whose value is of the order of the Fermi energy of the protons and neutrons when 
their mean separation is of the order of a pion Compton length, i.e. E nuc ~ (m 7r c) 2 /2m p . 
Above this density, 




which works out to be a few times 10 11 gm/cm 3 , the crust matter will still contain positively 
charged baryonic nuclei in a negatively charged Fermi sea of electrons, but there will now also 
be a third constituent consisting of freely moving neutrons outside the nuclei. 

The use of the term "crust" to describe the layers both above and below the critical 
value (||) is motivated by the consideration that the ionic nuclei will crystallise as a Coulomb 
lattice whose large scale behaviour will be that of an elastic solid as soon as the star has 
cooled sufficiently Except for a very thin outer surface layer with density below the white 
dwarf limit (H) that may remain liquid as a relatively shallow "ocean" , the rest of the crust is 
expected Q to have solidified by the time the neutron star has dropped below the MeV level, 
which will be reached within a few months of its formation. Due to the high conductivity of 
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the degenerate electrons the outer magnetic field will be firmly anchored in this crust, whose 
rotation rate is therefore what is is measured directly by pulsar frequency observations. 

For the purpose of explaining these observations (see Figure ||) the most interesting fea- 
ture of the crust is the presence of the interprenetrating neutron fluid in the inner crust, at 
densities ranging from the critical "drip" value @ all the way up to the nuclear value (||) that 
is reached at the base of the crust beyond which the ions dissolve. In the relatively low tem- 
peratures (below an MeV) that are relevant it is generally believed [J that the unconfined 
neutrons (like those within the nuclei) will form Cooper type pairs that will form a bosonic 
condensate. The interpenetrating neutron constituent is thereby endowed with the property 
of superfluidity, which enables it to flow freely past the metallic lattice (and the electron sea 
to which the lattice is electrically coupled) in the manner illustrated schematically in the fol- 
lowing diagram (using hyphens to indicate the negatively charged electrons, crossed circles to 
indicate the positively charged nuclei, and double arrows to indicate relatively moving Cooper 
pairs of neutrons): 

- e e - 
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The unconfined neutrons thus constitute a massive component that can rotate independently 
of the crust, thereby - as explained below - providing the most promising kind of mechanism 
for explaining the observed pulsar frequency glitches (see Figure [[]). 

At the base of the inner crust, at densities above the nuclear value given by (||), it is 
generally believed that the neutron fluid and ionic constituents merge to form a uniform 
fluid composed mainly of (superfluid) neutrons but with an independently moving (supercon- 



ducting) protonic constituent, and with the further complication [ 10 1 that instead of forming 
ordinary scalar Cooper type pairs the neutrons at this deeper level condense as pairs of spin 
1. At even higher densities, beyond that of ordinary nuclear matter, various more or less 
exotic possibilities have been suggested, but no firm concensus has yet emerged. For exam- 



ple Glendenning has predicted [11, 12] that there will be a hybrid zone in which negatively 
charged drops of quark matter will condense within the surrounding positively charged bary- 
onic liquid, and moreover that they will crystallise to form an ionic solid analogous to that of 
the crust. At even higher densities, as the maximum allowed by (||) is approached, one might 
expect that there would be an inner core where the drops merge to form another homogeneous 
superconducting superfluid zone, that (unlike the outer, baryonic core) would be constituted 
purely of quark matter, in which interesting new kinds of superfluidity and superconductiv- 
ity @, m, u , n c °u id 

occur. Like the somewhat better understood inner crust and outer 
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Figure 1: Qualitative sketch of a typical observational plot of pulsar angular velocity £1 against 
time t. The long down - pointing arrow indicates the negative change AO during a period 
of steady slow down. The short up - pointing arrow indicates (on an exagerated scale) the 
positive jump 5Q during a glitch (consisting of a sharp discontinuity followed by a transient 
readjustment). 

core regions, these very high density inner regions may also be relevant to the interpretation 
of pulsar frequency observations such as are illustrated in Figure [|. 

The overall situation is not just that the global structure and behaviour of a neutron 
star is rather complicated, but furthermore that (as a fortiori for strange stars if they exist) 
many important aspects remain so unclear that - except in the crust region (or for very 
low mass neutron stars) for which a reasonable degree of concensus already prevails - it is 
hardly worthwhile yet even to start the detailed numerical calculations that will be needed 
later on. Before a convincingly realistic neutron star model can be developed even as a 
rough approximation many underlying physical issues will need to be dealt with, of which the 
most basic are concerned with the qualitative nature of matter at the supernuclear densities 
attained in the cores of all neutron stars above or near the precise Chandrasekhar limit 
value, M ~ \/2M a (where M Q is the solar mass) though not for very low mass neutron star 
configurations (which are at least of academic interest, even though it is hard to see how they 
could be created in nature. 

From the point of view of the interpretation of observational data, many less fundamental 
but technically non-trivial issues need to be be clarified. Among the other accessory issues 
(concerning matter at less extreme densities) that also need to be dealt with, the one with 
which the present article will be primarily concerned is that of the consequences of the pre- 
dicted superfluidity. The final sections will concentrate on the results of recent progress on 
the development of an appropriately relativistic treatment as an improvement (in view of 
immediate coherence, as well as the long term objecive of precision) over the non-relativistic 
treatment that has until now been mainly used, not just for superfluidity but also for many 
other relevant phenomena. However before going into the technical aspects of the relativistic 
treatment, the first part of this article will provide a brief survey of the reasons why the 
phenomenon of superfluidity is particularly important for relating theoretical understanding 
of the inner structure of the neutron star to the available observational data, of which the 
most richly informative part (see Figure |TJ) comes from pulsar timing. 
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1.2 Minimal two component rotating star models 

As emphasised above, for many of the most important questions about the global structure 
of neutron stars, no quantitative agreement is available or to be expected for a long time 
yet. There are however several essential qualitative features on which practically all neutron 
star theorists do seem to agree already. In relation to the pulsar timing observations, the 
most important of these agreed features J?]] is the presence of a rigidly corotating structure 
that may or may not include part of the core but certainly includes a solid outer crust to 
which is anchored the magnetic field configuration that gives rise to the observed radio and 
other emission. The next most important feature, common to all viable theoretical scenarios 
albeit of a more subtle nature, is the presence of some (maybe many) effectively superfluid 
(or superconducting) zones that can rotate independently ||] of the rigidly rotating crust 
structure whose angular velocity, Q c say, is presumed to be the same as the f2 that is directly 
observed. 

The mimimal agreement about these two essential features is what underlies the long- 
standing, widespread, and enduring popularity of a corresponding kind minimally complicated 
rotating neutron star model, involving just two independently rotating parts: a "corotating 
crust" part with (directly observed) angular velocity f2 c and a "superfluid neutron" part with 
a possibly different angular velocity f2 n (representing the average of what in a more detailed 
treatment would be a spacially variable angular velocity distribution). 

The basic postulate of such a minimal model is that the total angular momentum J of the 
star is the sum of decoupled parts, 

J = Jc + Jn , (6) 

with 

J c — Iq^c ; Jn — IjaQc j (7) 

where I c and I n are separate moments of inertia that are supposed to remain constant during 
a process of continued variation governed by an external torque r ex . (A more sophisticated 
variant of this model would allow for cross coupling, whereby J c is affected by Q n and vice 
versa: a small cross coupling of this kind would inevitably be present [17] in an exactly 



relativistic description, and a possibly more important cross coupling effect is to be expected 
from the effect - to be discussed below - of superfluid momentum "entrainment" , whose 
likely relevance in neutron star matter was originally pointed out in the context of proton 
super conduct vity |l9]|.) 

The external couple T ex represents the effect of the magnetic field anchored in the rotating 
crust, which, in view of the high conductivity of the crust is generally assumed to remain 
constant over timescales long compared with those (at most a few years, since the oldest 
pulsar observations go back only to 1968) of the observed fluctations. The effect of this 
steady couple is of course to cause a total angular momentum loss rate given by 

J = r cx < o . (8) 

If the angular velocities were locked together, Q n = Q c = Q, this would give f2 = r ex /I with 

I = I c + In, (9) 

and this relation will in any case be true for the long term average, {Vt) i.e. since (for reasons 
to be discussed below) the separate angular velocities can never get too far appart, the long 
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term slowdown of the observed pulsar frequency allows the torque involved to be estimated 



as 



j<n>. (io) 



in which, for neutron stars with mass M ~ 3M /2 (which, consistently with the theoretical 
prediction (|]), is what has been found |2Q] for the few cases in which the mass is reliably 
measurable) , the total moment of inertia (unlike the distinct parts I n and I c for which different 
theoreticians have rather diverse ideas in various cases) can be evaluated in a generally agreed 



manner, which leads [21] to estimates of about 10 M Km within a factor of order unity 
(whereas the uncertainty range would be much larger for a neutron star nearer the upper 
mass limit). 

The idea of the two component model is that as well as supporting the effect of the external 
torque r ex , the crust component exerts an internal torque Ti n on the "superfluid neutron" 
component, which therefore obeys an evolution equation of the form 

Jn = Tin , (11) 

while, in order to be consistent with (||) the crust component must obey 

Jc — -Tex Tin j i\^J 

in which, unlike the external torque r ex , the internal torque is not constant but proportional 
to the angular velocity difference 

Tin = -^r^-, UJ = tl n -Q c , (13) 

1 T 

where r is a damping timescale whose estimation will be discussed below. The chosen nor- 
malisation of this timescale is such that, according to the preceeding equations, the angular 
velocity difference will satisfy a differential equation of the simple form 

UJ + - = -y(n), (14) 

in which the right hand side is a constant that fixes the saturation limit 

u —* uj s , uj s = — — (tl) t > , (15) 

to which the difference lv will tend in the long run, unless this continuous evolution process is 
interrupted by a "glitch". So long as no such interruption has occurred, the angular velocity 
difference will be given as a function of the time t by an expression of the form 

oj = lu s + (oi - w s )e~ t/r , (16) 

where loq is a constant of integration interpretable as the value of to when the clock time t 
was set to zero. 
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1.3 The problem of accounting for glitches 

Some of the strongest observational evidence in favour of the theoretical picture epitomised 
by the kind of highly simplified 2 component neutron star model described in the previous 
section is provided by a phenomenon in which the continuous evolution described by this 
model is subject to a temporary model break down. The phenomenon in question is what is 
known as a glitch (see Figure [[]) during which (with a rise time too short to be measured, at 
most a few hours and probably much less) the observed angular frequency 0, of an isolated 
pulsar undergoes a positive discontinuity, > 0, that partially cancels the loss AQ during 
the preceeding period of continuous slowdown. 

Since there is no imaginable way the external torque could suddenly become very large (nor 
any observational evidence that the associated pulsar emission process changes significantly 
at all during a glitch) there can be no corresponding discontinuity in the angular momentum. 
This means that, if we want to use a model only involving a single component, it will be 
necessary to take account of variation of the moment of inertia, whose total / will evidently 
undergo a negative variation 51 given, for a glitch of amplitude 50, by 

AO 

5J = 81 = -I (17) 

The earliest theory designed to account for this phenomenon, as first observed in the Crab 
and Vela pulsars, was based on the first of what were presented in the preceeding section as 
generally agreed features of neutron stars, namely the presence of a solid crust structure, but 
on the basis of a single component model taking no account of the second generally agreed 



feature (namely the possibility of independently rotating parts). The idea [22, g, |23| was 
that the rigidity of the solid crust would tend to prevent the decrease in moment of inertia 
that would necessarily accompany the loss of angular momentum in a purely fluid star. In a 
simple rotating fluid star model, the oblateness due to centrifugal force would give rise to a 
variable moment of inertia that would be expressible for low values of the angular velocity by 
an expression of the quadratic form 

n 2 



J - J °( 1 + ns)' ( 18 ) 



where Jo is the value of the moment of inertia in the non-rotating spherical limit and f2* is 
a constant specifying the relatively high value of the angular velocity (which will be given 
in terms of the mean density subject to @ by the rough order of magnitude estimate ~ 
G/9*) that would be necessary for deviations from spherical symmetry to be of order unity. 
Whatever may have happened immediately after the birth of the neutron star, no such rapid 
rotation still occurs in any of the (at least centuries old) pulsars that are actually observed 
today, for which the condition f2 2 <C fi* is always satisfied. 

For a simple perfect fluid star model, the effect of the external torque (||) during an 
extended time interval At would be to cause an angular momentum loss, AJ ~ r ex At < 
that would be accompanied by a corresponding angular velocity variation AS7 < 0, which 
according to (|l8|) would entail a decrease in moment of inertia given by 

n 2 An 

AI ^ 2I w^r <0 - (19) 



8 



Due to the solidity of the crust, which tends to preserve the more highly elliptic initial config- 
uration, the actual change in the moment of inertia will fall short of what is predicted by this 
formula, but at some stage the strain will build up to the point at which the solid structure 
will break down. The implication is that there will then be a "crustquake" , in which the solid 
structure suddenly changes towards what the perfect fluid structure would have been, thereby 
changing the moment of inertia by an amount that will be at most of the order of the upper 
limit given by 

SI < A/, (20) 

where AI is what is given by (|l9| ), and that will be considerably less than this if the crust 
rigidity is low. According to (|17| ) the corresponding positive angular velocity discontinuity 
<5$7 associated with the continuous negative angular velocity change Af2 since the preceeding 
glitch will be subject to the limit 

6Q<-2^An, (21) 

The preceeding formula provides an order of magnitude limit that must be satisfied by a 
rather large margin if the rigidity is low but that is entirely consistent with what is observed 
in the case of the Crab pulsar, for which typical glitches are characterised by 5tl ~ 10 _8 fL 
However almost immediately after it was first proposed, it began to be recognised @ that 
this rather obvious mechanism would not be sufficient to account for the much larger glitches 
that are frequently observed in cases such as that of the Vela pulsar, for which typical glitch 
amplitudes are characterised by 50, ~ 10 -6 fi ~ — 10 _2 Afl. 



Since it was first suggested by pioneers such as Anderson and Itoh [24], the generally 
accepted way of getting round this limitation - namely that the likely changes of the moments 
of inertia will be far too small to account for frequent giant (Vela type) glitches - is to drop 
the single component description in favor of the two component description in which glitches 
can be accounted for even if (as assumed in its simplest version) the relevant moments of 
inertia undergo no significant change at all. The essential point is that the consideration that 
the very short glitch duration excludes any significant jump in the total angular momentum 
does not rule out the possibility of impulsive transfer of angular momentum between the two 
components provided they balance out: 

5 J = SJ C = -5J n , (22) 

so that 

Sn Q = -ksn c . (23) 

The idea is that between the glitches the weak coupling mechanism described by (|TJ) allows 
the slowdown of the "neutron superfluid" angular velocity f2 n to lag behind that of the crust 
component which is what is presumed to be actually observed, f2 = f2 c , so that during the 
preceeding period At the angular velocity difference lo will be positive. It is generally supposed 
that (for diverse reasons to be discussed below) this angular velocity difference gives rise to 
stresses that are partially relaxed in the glitch process, whose onset occurs when the difference 
lo reaches a critical glitch inducing value ui g say that, in order to be attainable must be less 



than the limit u; s given by (15) - a condition that would fail if the relaxation timescale r were 
too short. 
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Leaving aside cases for which u g > lo s (whose evolution will be of the glitch free kind 
recently investigated |^] by Sedakian and Cordes) as well as the marginal case in which 
Lj g u> s , i.e. subject to the proviso that there is a safe margin uj g <C u s , the evolution 
equation ( |l6| ) will be replacable by the linear relation 

^ = ^ + L, (24 ) 

in which each of the terms is small compared with unity. Successive glitches bring about 
negative adjustments 800 that are needed to cancel out the cumulative effect of the positive 
variations A that develop during the duration of the interglitch periods governed by (^4|), so 
that on average they cannot deviate too much from the order of magnitude estimate given 
simply by 

5uj^-Auj, (25) 



in which, by (15) and (24), the deviation built up during an interglitch interval of duration 



At will be given simply by 

Au = - — (n)At. (26) 

Using (|23| ) to eliminate the unobservable jump <5S7 n from the difference 5oj = 5£l n — 5fl c , the 
magnitude of the observable jump 5£l c can be estimated by (p5|) as 

<50 c ^h-Acu. (27) 

Since the observable interglitch frequency variation will be given roughly by AQ C ~ (17) At 



one sees from (26) and (27) that it provides a corresponding estimate 



5n c ^-—An c , (28) 
for the observable frequency jump during a glitch. The presumption that S7 C is identifiable 



with the O that is observed allows us to compare this with the previous upper limit (21) 
that was obtained for the single component model with variable moment of inertia. It can 
be seen that the difference is simply that the small factor (fJ/fi*) 2 in the upper limit for the 
single component model is replaced, in the two component fixed moment of inertia model, by 
the ratio I n /I c whose value is highly uncertain (in view of our lack of firm knowledge about 
what goes on in the core of the neutron star) but can plausibly be supposed to be of the 
order of unity, which is what is needed to account for the frequent very large glitches (with 
<5$7 « — 10~ 2 A$7) that are observed in examples such as Vela. 



1.4 The question of pinning and the damping timescale 

The foregoing estimate (^8|) is not sensitive to the particular value of the damping timescale r 
except that it is assumed to be large compared with the interglitch period At which is usually 
several months or more. This requirement might at first sight appear to be incompatible with 
observations of post glitch relaxation, in which shorter timescales of only a few weeks have 
been shown to be involved. Such discrepancies are however to be expected on the basis of 
our general qualitative understanding |26|] of what is involved. The strong density gradients 
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in the star imply the existence of many different zones in which differential rotation with a 
wide range of damping timescales can occur. Our simplified two component description of 
glitches depends on taking the part with moment of inertia I n to correspond to a substructure 
for which the relevant damping timescale timescale r is very large. A formally similar two 
constituent model might also be used for describing post glitch relaxation with much shorter 
timescales, but for such an application the substructure with moment of inertia I n would have 
to be reinterpreted as corresponding to some other part of the star. Of course if we wanted 
to describe both the glitches and the postglitch relaxation, in a single coherent framework, 
we would need to amalgamate the separate two component models so as to obtain a more 
elaborate model (such as has recently been used [27] for the analysis of precession) with 
three or more independent components (and with not just a single damping timescale but 
an antisymmetric matrix of mutual damping coefficients). Although the construction of such 
composite models is straightforward in principle, most authors have so far (quite reasonably) 
preferred to concentrate on particular aspects for which a less intricate description is adequate. 

Even for applications, such as the glitch model of the previous section, for which a two 
component description is adequate as a lowest order approximation, the estimation of the 
relevant damping timescale remains a subject of great uncertainty. The situation has however 
been clearer since the general question of quasi stationary equilibrium in a rotating superfluid 
was systematically addressed in the context of neutron stars by Alpar and Sauls [28], Bildsten 



and Epstein |2{|] and the Sedrakians [ 30 1 , who drew attention to the consideration that long 
damping timescales can arise not just from weak but also from strong coupling. These authors 
considered the basic general problem of a two constituent superfluid model of the simplest 
kind in a local configuration of differential rotation about a fixed axis characterised by a 
unit 3 vector, V say in the neighbourhood of a position determined (in a Newtonian flat space 
description) with respect to a central rotation axis by an orthogonal radius vector r. One of the 
constituents is the "normal" (and therefore in a n equilibrium state) rigidly rotating) crust 
constituent characterised by a (uniform) angular velocity Cl c and a corresponding velocity 
vector given as the cross product v c = fl c V x r. The large scale averaged velocity of the 
superfluid constituent - which for our purpose is to be thought of as constituted of neutrons - 
is characterised in terms of a perhaps radially variable angular velocity Sl n by a similar formula 
v n = f2 n v x f. However in the latter case it is to be born in mind that that on a microscopic 
scale the superfluid fluid is irrotational except on quantised vortex lines round which the 
integral of the relevant superfluid particle momentum mv n is given by the Planck constant, 
i.e. 2irh in Dirac's notation, so that the corresponding velocity circulation is K = 2irh/m - 
while, in view of the bosonic pairing, the relevant mass scale in the application with which 
we are concerned is twice that of the neutron, i.e. m = 2m n 

Since the averaged velocity circulation per unit area orthogonal to the rotation direction 
v will simply be 2Q n , it follows that the corresponding surface number density, a say, of 
quantised vortices will be given by 

Cirri , . 

a = — . (29) 
nh 

Although the superfluid motion is non dissipative, and so has no direct interaction with the 
"normal" crust material, the vortex cores (which are defects where the superflidity breaks 
down) will in general be subject to a drag force -Fd say per unit length, exerted by the 
bacground in the direction of relative motion. Using the notation v v for the velocity of motion 
of the vortex lines orthogonally to the rotationaxis V the dissipative drag force exerted by the 
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background will be given by a formula of the form 



F d = 7] T (v c - v v ) , (30) 

in which r] T is a positive resistive drag coefficient, whose quantitative evaluation is a subject 
of much uncertainty - not just in the core, but even in the crust, where it is very sensitive to 
temperature and other quantities such as superfluid pairing correlation lengths that are 
rather difficult to estimate j p2"[ . 

This drag will not be the only force acting on the vortex line, which will also be subject to 
the Magnus effect. According to the well known formula J| of Joukowski (or, in an alternative 
transliteration, Zhukovskii) this gives rise to a lift force F\ per unit length that is proportional 
to the product of the particle number density and the corresponding momentum circulation: 

F\ = 2irhn(v n — v v ) x u, (31) 

where n is the relevant particle number density. In the neutron superfluid application with 
which we are concerned the effect of the bosonic pairing must be taken into account, which 
means that the relevant number density is only half that of the neutrons themselves, i.e. 
n = n n /2. 

On the microscopically very long timescales characterising the relevant applications the 
vortex lines can be treated as effectively massless, which means that the evolution of the 
system will be determined simply by the condition that the total force on a vortex line must 
cancel out, 

F d + F l = 0. (32) 

To solve this, it is convenient to decompose the velocity v v of the vortex lines into a (small) 
radially outward directed part, with magnitude r, and a (larger, but for our purpose less 
important) remainder, directed parallel to the fluid flow vectors, in the form 

j. 

v v = - r + U v v x r, (33) 
r 

where fi v is interpretable as the angular velocity of the vortex lattice. It is also convenient to 
introduce a dimensionless resistivity coefficient defined by 

"* (34) 



27rhn ' 



which is what in the jargon of aero engineering would be called the drag to lift ratio (what, in 
that context, one seeks to minimise by cunning aerofoil design). The solution of (^) is thereby 
expressible as the condition that the vortex line angular velocity is intermediate between those 
of the crust and supefluid constituents, with value given by 

O -4- r 2 Q 

V v = ^ + Cr " c , (35) 
1 + cf 

while the radially outward "creep" component of the velocity of the vortex lines will be given 
by 

r = (n n - n c ) . (36) 

1 + ci 
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This last equation is particularly important because it determines the rate of change of 
the vortex line surface density a: as the comoving radius of the vortex distribution increases, 
the surface density will evidently undergo a corresponding decrease given by the relation 
a jo = —2r/r. Since by ( |29| ) this surface density is proportional to the superfluid angular 
velocity, we deduce that the rate of variation of the latter will be given by 

h n _ 2c r u; 

where u is the angular velocity difference as introduced in the relation (|i~3|), from which, by 
comparison with (37), the corresponding value of the damping timescale r in which we are 



interested in, can be read out as 



fifc(* + £)- (38) 



T : 2m 

As well as showing that the timescale is subject to a lower limit (that might have been 
guessed on dimensional grounds) given by r > I c /IQ n and attained for c r = 1, a noteworthy 
feature of this result [g8|, ^9], |30|, is the dual symmetry between the roles of the drag to lift 
ratio c r and of its inverse, the lift to drag ratio c~ l . The decay timescale r becomes infinitely 
large not just in the drag free limit for which rj8r and hence c r become arbitrarily small - so 
that by (|35|) the vortices are dragged allong with the superfluid - but also in the opposite 
"pinned" limit of very large c, for which the force on the vortices is strong enough to lock 
them to the crust material. 

The estimation of the actual values of the c r , in the various zones of interest, has been the 
subject of much work, but the subject is difficult and many of the results are still inconclusive 
or controversial. Following the recognition p^] that magnetic coupling forces between the crust 
and the superconducting proton neutron superfluid zone are more important than had been 
previously supposed, and strong enough to lock this core region to the crust on timescales short 



compared with those relevant to glitch observations, it was suggested [33] that even the new 
coupling force values were underestimated, so much so that c r would become large compared 
with unity, with the implication that r could be very long after all, just as had been supposed 
in the early years when c r had been supposed to be small. Even in the qualitatively more 
familiar crust regime the situation is still unclear, partly because of effects of temperature 
dependence: work of Jones |3^, |J] suggests that pinning may be much less effective than had 
been previously supposed so that instead of being high c r would be very low there. 

1.5 The long term crustal drift phenomenon 

The question of the effectiveness of vortex pinning to the crust comonent leads on to the 
related issue of what is actually responsible for the stress whose release, when a critical 
value is exceeded, is supposed to provide the glitch mechanism in the two constituent scenario 
described in the preceeding section. All the early versions of such a two constituent mechanism 
assumed that the relevant stresses would be due to vortex pinning. Like their more recent 
variants, the various early versions were classifiable in two distinct categories. In the first 
category p3, 35] it was supposed that the discontinuous breakdown would occur when some 



maximum static pinning force value was exceeded: the sudden (rather than "creeping" ) nature 
of the breakdown was accounted for, in a recent version |36|], as being due to a thermal 
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instability resulting from the temperature sensitivity of c r , while another new suggestion [37] 
is that the relevant slippage occurs at the locus where the ions dissolve at the base of the 
crust. In the second category [^4|, it was suggested that such a maximum pinning 

force value might never be reached because the elastic solid structure would breakdown first 
in a crustquake (of the kind required in the single constituent moment of inertia changing 
mechanism that may account for cases such as that of the Crab). 

In all these various versions, the necessary transfer of angular momentum from the relevant 
independently rotating layers with moment of inertia I n to the crust component with moment 
of inertia I c is mainly attributable to the torque exerted by the pinning forces. However it 
has recently been pointed out [E]] that there is an alternative possibility (effectively a new 
variant within the second category) whereby the necessary angular momentum transfer may 
be acheived convectively - by a transfer of matter (removal from the crust of matter with low 
angular momentum, and its replacement by matter with higher angular angular momentum) 
which can occur even if torque forces are entirely absent, i.e. in the small c r limit. In this new 
kind of scenario, the stress ultimately responsible (when a critical level has been exceeded) 
for the discontinuous transfer is attributable to a centrifugal buoyancy deficit in the relatively 
slowed down crust component. 

In the case of the earlier pinning driven mechanism, it was pointed out by Ruderman |39[| 
that if the glitches were due to breakdown of the solid structure (rather than discontinuous 
vortex slippage) then the long term effect of many glitches would be analogous to that of 



terrestrial contentinental drift. It would give rise to a pattern of convective circulation [41] 
involving "transfusion" of matter from the crust constituent to the underlying neutron super- 
fluid constituent in a "subduction" region near the equator, at colatitude 6 = ir/2, and the 
other way round near the poles at colatitude 8 = 0. The corresponding long term average 
rate {9) of angular drift of a crust plate at the surface, which for the spin down of an iso- 
lated pulsar would be directed away from the pole towards the equator (see Figure ^), was 
estimated by Ruderman on the assumption that it would correspond to an outward velocity 
of the same order of magnitude as the mean cylindrical expansion rate, (r) ~ —(&)/2a of the 
vortex distribution, whose surface number density a is given in terms of the angular velocity 
by the proportionality relation (p9[). This reasonning |39] provided a formula of the form 



<*— (39) 

which implies that the timescale for complete turnover of the crust material is of the same 
order as the spin down lifetime of the pulsar, during which, as Ruderman pointed out, the 
magnetic dipole would be able to be dragged most of the way from the rotation axis to the 
equator. This would result in a net increase (ocurring discontinuously at the glitches) of the 
pulsar radiation rate and thus of the magnitude of the spin down rate Cl. Another reason for 
an increase of the spin down rate would be the decrease in oblateness according to (|l8|), but 
this would evidently be much less important. 

Unlike the original single constituent mechanism |2^, ||, ^3], [7j based on the loss of moment 
of inertia due to decrease in oblateness, and unlike the versions [35, 36, [37]] of the two- 
constituent theory that attribute the glitches to discontinuous vortex slippage, but like the 
Ruderman version p8| , |39| ] (that applies when the pinning is too strong to be broken) the 
newly proposed two-constituent mechanism ^(| (that applies when the pinning is too weak 
to be effective) will also entail a substantial rate of long term drift of plates of crust material. 
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Figure 2: Qualitative sketch indicating direction of force expected (c.f Ruderman 1991) to act 
on (magnetically slowed down) crust due to vortex pinning mechanism, if it is effective, when 
the (interpenetrating) neutron superfluid retains a higher rotation rate. (Vertical shading 
indicates the alignement of the vortices in the region occupied by neutron superfluid, which 
is not confined to the core but interpenetrates the greater part of the solid cust as well.) 



However this centrifugal buoyancy deficit mechanism differs from Ruderman's pinning driven 
mechanism in a manner that may be experimentally observable, since it is expected to produce 
plate drift in just the opposite direction, meaning that of decreasing colatitude 9 for an isolated 
spinning down pulsar (see Figure |3j) , entailing transfusion of matter into the crust constituent 
near the equator, and out of it nearer the poles where 9 is small. In this mechanism (unlike 
Ruderman's) the angular momentum of a crust plate will not be significantly changed when 
its colatitude undergoes a displacement 59 during a glitch, so its change of rotation frequency 
can be estimated as being given roughly by 5Q/20, ~ —59 where 5£l is the glitch amplitude 
that is actually observed, and that partially cancels the preceeding interglitch variation AQ. 
The change observed in the long run is the sum over the glitches of the combination Af2 + 5$l, 
which will be the same as the sum of the hidden changes 5 n (since the interglitch variation 
A n of the relevant superfluid part is assumed to be negligible). Since 5Q n = —(I c /I n )5Q, by 
(|23|), it can be seen to follow that the long run average of the angular drift rate will be given 

by 



which has the opposite sign to what is given by the Ruderman formula fl39|), but can be 
comparable in magnitude since I n can be comparable with I c . Indeed, in a case for which 
the moment of inertia I n of the hidden part is large compared with the crust contribution, 
the magnitude given by the new formula (40) would be correspondingly larger than in the 
previous case, with the implication that the crust material would be entirely recycled several 
times during the spin down lifetime of the star, while this lifetime itself would presumably be 
considerably prolonged because the magnetic dipole axis would be dragged towards the pole, 
thereby decreasing the pulsar radiation rate. 

For the purpose of observational discrimination between cases involving strong [38, 3£], 



moderate [35, 36 1, or very weak pOj| , coupling the relevant directly measurable parameter is 
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Figure 3: Qualitative sketch (using same shading conventions as before) indicating direction of 
force expected to act on (magnetically slowed down) crust, even if vortex pinning is ineffective, 
due to centrifugal buoyancy mechanism when the (interpenetrating) neutron superfluid retains 
a higher rotation rate. 



what is known [42] as the (long term) braking index, 



n=m, (41) 

(ft) 2 ' V ' 

and more particularly the braking deficit 

e = 3 - n (42) 



between the value that is observed and the value, n = 3, that is predicted |43j| for a simple, 
rigid, non aligned magnetic dipole model, and also for more sophisticated pulsar emission 



models including allowance 44] for outflow of charged particles. If it is assumed that particle 
outflow and changes of moment of inertia can be neglected, then according to the simple 
dipole model ][E|] the relative spin down rate Cl/Q is just proportional to (fisina) 2 where a 
is the dipole misalignment angle, i.e. the colatitude of the magnetic pole, so for this case the 
difference (42) can be immediately evaluated as 



2!! <"> (43) 



tana (Cl) 



When it gets near the extreme polar or equatorial values, a ~ 0, or a ~ tt/2 , the evolution of 
the misalignment angle will of course have to come to a halt, (a) — 0, but in the intermediate 
range, i.e. for tana ~ 1 one would expect the misalignment angle to move with the crustal 
drift, (a) ~ (9). Subject to this assumption, the Ruderman formula ( j39[ ) for the strong 
pinning model leads to the positive estimate 

6 ~ cot ol , (44) 

while the formula (^) for the model with negligible pinning gives the negative estimate 

e — -p-cot a . (45) 
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The effects of variation of the moment of inertia |49] should not substantially effect the 
validity of these estimates except for new born pulsars with extremely rapid rotation, but 
electromagnetic effects of various kinds f46| , 47] (including the obvious possibility of magnetic 
field decay) are more likely to give significant, typically positive, contributons to e. (In glitch 
free scenarios, Sedrakian and Cordes [25] have pointed out that differential rotation may bring 
e down to negative values for periods of limited duration, but this sort of effect can be expected 
to cancel out in a long term average over many gliches). 

If these estimates are indeed applicable, then in the case of the Vela pulsar the observed 
value pS[] , namely e ~ 3/2, can be plausibly construed as evidence favouring the Ruderman 
model, with a moderate misalignement angle of the order of 40 degrees. A less clear cut 
case example is that of the Crab, for which the observed value J48|], namely e ~ 1/2, is also 
positive but considerably smaller, which suggests that this may be another instance to which 
Ruderman model applies, though with a relatively high misalignment angle. However in view 
of the above mentionned likelihood [{461 , (4^ 1 of other of other positive contributions to e, this 
evidence is too inconclusive to exclude the possibility that the Crab glitches may, after all, 



be attributable a slippage mechanism |35, |37|] of the first category, or even to the original 



simple oblateness mechanism subject to (21). 

What transpires from all the work that has been rather rapidly surveyed in the proceed- 
ing sections is that the available theory of the internal structure of neutron stars seems to 
provide all the elements needed to account for the accumulated pulsar frequency data within 
the framework of scenarios in which superfluidity and differential rotation commonly play an 
essential role. However we are still a long way short of reaching any concensus about the 
detailed modelling of the many different kinds of behaviour observed in particular cases such 
as the Crab and Vela pulsars. Before any definitive understanding can be reached it will 
be necessary to carry out much more work on the technicalities of basic physical processes, 
particularly those involving electromagnetic effects, which were barely mentioned in the pro- 
ceeding overview, but that are extremely important for the detailed estimation of important 
quantities such as the drag to lift ratio, c r . 

The remainder of this article will be concerned with just one of the many technical 
problems that needs to be dealt with before a satisfactorily complete understanding can 
be achieved. This is the problem of developing an appropriately relativistic treatment of the 
superfluidity that has been seen to play such an essential role in accounting for cases such as 
that of Vela, and even the less extreme case of the Crab. 



2 Essentials of relativistic superfluid mechanics 
2.1 Motivation and background. 

In calculations of global quantities such as the mass and radius of a neutron star with a 
given baryon number, it has been known since before the earliest pulsar observations in 1968 
that a fully relativistic treatment is indispensible for even a minimally acceptable level of 
accuracy. It is fortunate that quantities such as this j|] can be obtained within the framework 
of an exactly spherical perfect fluid description for which a fully relativistic treatment is easily 
applicable and has always been used. This contrasts with what has been done about secondary 
effects such as precession [50, 27, 51], involving mechanisms such as elasticity (for which a 



relativistic treatment has long been available |52|, |53J but is difficult to apply) and superfluidity 
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(for which the relevant macroscopic treatment |^4|, |5(| is relatively new). Even in the 
relatively tractable context of stationary axisymmetric problems, or in contexts involving the 
(intrinsically relativistic) effect of gravitational radiation, superfluidity and superconductivity 
have nearly always been dealt with using a non-relativistic Newtonian, even in relatively recent 
work © H H 

There are two essentially different reasons why it is worthwhile to try to do better. One 
is of course that a fully relativistic treatment should in principle be more accurate, and 
will no doubt become necessary for this purpose sooner or later. However in the short run 
this is not always what is most important, since the errors inherent in the use of a purely 
Newtonian treatment range typically from a few per cent to a few tens of percent which is 
not very significant compared with the order of unity (or worse) uncertainties about many 
of the physical quantities involved. The other kind of reason (which some readers may find 
surprising) is that for many purposes the use of a relativistic treatment is actually easier. 
In many cases the advantage of a relativistic treatment is due to the fact that the Lorentz 
group is in the technical sense semi-simple, whereas the Galilei group unfortunately is not. 
However, whether or not it is intrinsically simpler, the use of a relativistic treatment will 
usually be more convenient in practice whenever one wishes to use the commonly appropriate 
strategy |2l|, |6^, of working in terms of perturbations of the available spherically symmetric 
perfect fluid neutron star models. This is just because (as remarked above) all the best of these 
models (the only ones that are commonly taken seriously) are already formulated exclusively in 
a (general) relativisic framework. The same consideration applies to the perturbations of the 



relativistic axisymmetric rapidly rotating star models that have recently |62|, 63, 64, 65, |6q ] 
been a subject of rapid development, particularly in relation to the question of bar mode 
instabilities that may be significant as a source of gravitational radiation. 

The treatment provided here will be limited to the case of scalar (spin 0) models such 
as are appropriate for the experimentally familiar example of helium - 4 (though not |67j| 
helium - 3) and also for the mesoscopic description of the neutron fluid that (as discussed 
above) is predicted to interpenetrate the ions in the lower crust of a neutron star, and that is 
believed |68| to condense as a superfluid in which the relevent bosons are scalar Cooper type 
pairs of neutrons. For the mesoscopic (intervortex) treatment of the mixed proton neutron 
superfluid below the crust a scalar description has commonly been employed |69| in a 
Newtonian treatment, and an analogous relativistic description [ff(], has recently been 
made available. However for an exact description of such a mixed proton neutron superfluid, 
in which it is predicted (6^] that the neutrons pair up as bosons of spin 1, it would be necessary 
to use a more elaborate treatment that has yet to be developed). 

What is actually needed for the analysis of large scale effects (such as were considered in 
the preceeding sections) is not just a mesoscopic treatment of the superfluid on scales small 
compared with the spacing between interpenetrating ionic nuclei and the vortices where the 
irrotationality condition breaks down, but a macroscopic average over much larger scales. An 
appropriate macroscopic theory of the kind that is needed has recently been developed |5(| and 
is described in the final subsections of this article. The treatment presented here differs from 
the relativistic analogue [^4|, 55] of the earlier non relativistic description of the averaged 
effect of vortices in neglecting the small anisotropy due to their effective tension, but instead it 
includes allowance for what in the long run is likely to be a more important effect, namely the 
"transfusion" of matter (for the reasons discussed in the previous subsections) between the 
superfluid constituent and the normal background, which is to be interpreted as representing 
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not just thermal excitations (as in ordinary liquid helium - 4) but the entire crust component. 
A more elaborate treatment would include allowance for anisotropy of the crust constituent 
which, as noted above, will be cold enough (except very near the surface) to behave as an 
elastic solid: the way to do this has been indicated elsewhere [73], but such a mixed fluid solid 
description has not yet been developed in detail, and will not be dealt with in the introductory 
treatment provided here. 

As a preliminary to the construction of model [56| that is actually needed for the macro- 
scopic treatment of neutron star matter, this presentation starts by recapitulating the long 
well known essentials of the relativistic version of the single constituent kind of superfluid 
model that is appropriate for the description of Helium 4 at zero temperature, and of of the 
more recently developed generalisation [74, 76 1 to a two constituent model (of the kind 
whose non-relativistic analogue was originally developped by Landau) in which the second 
constituent represents a gas of phonon excitations. 



2.2 Single constituent perfect fluid models 

Before getting into the specific technicalities of superfluidity, it is worthwhile to start by 
recapitulating the relevant properties of ordinary barotropic fluid models, which includes the 
category of single constituent (scalar) superfluid models (representing the zero temperature 
limit of Landau type 2 constituent models) as the special case in which the vorticity is zero. 
The vorticity, in this context, is to be interpreted as the meaning the exterior derivative of 
the relevant momentum covector which will be formally defined below, so the vanishing of 
the vorticity is the condition for this momentum to be the gradient of a scalar potential, 
which in the superfluid case is to be understood to be proportional to the phase angle of an 
underlying bosonic quantum condensate. The qualification that this (zero temperature limit) 
model is barotropic simply means that there is only one independent state function such as 
the conserved (e.g. baryon) number density n or the mass density p (which are proportional 
in Newtonian theory but non-linearly related in relativistic theory) on which all the other 
state functions, such as the pressure P are dependent. The equation of state giving P as a 
function of p will also determine a corresponding speed c : say, of ordinary "first" sound, that 
will be given by the familiar formula 

c 2 = dP/dp , (46) 

and that must be subluminal, c 2 < c 2 (where c is the speed of light) in order for the usual 
causality requirement to be respected. 

Before proceeding it is desirable to recall the essential elements of the relativistic kinemat- 
ics and dynamics that will be required. This is particularly necessary in view of the regretable 
tradition in non-relativistic fluid theory - and particularly in non relativistic superfluid the- 
ory - of obscuring the essential distinction between velocity (which formally belongs in a 
tangent bundle) and momentum (which formally belongs in a cotangent bundle) despite the 
fact that the distinction is generally respected in other branches of non-relativistic condensed 
matter theory, such as solid state physics, where the possibility of non-alignment between 
the 3- velocity v a , and the effective 3-momentum p a of an electron travelling in a metallic 
lattice is well known. Failure to distinguish between contravariant entities (with "upstairs" 
indices) such as the velocity v a and covariant entities (with "downstairs" indices) such as the 
momentum p a is something that one can get away with in a non-relativistic treatment only 
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at a price that includes restriction to strictly Cartesian (rather than for example cylindrical 
or comoving) coordinates. 

In a relativistic treatment, even using coordinates <-> {t, x a } of Minkowski type, with a 
flat spacetime metric g^ v whose components are of the fixed standard form diag{— c 2 , 1, 1, 1}, 
the necessity of distinguishing between raised and lowered indices is inescapable. Thus for a 
trajectory parametrised by proper time r, the correponding unit tangent vector 

u» = ^f- (47) 
(it 

is automatically, by construction, a contravariant vector: its space components, u a = ^v a 
with 7 = (1 — f 2 /c 2 ) -1 / 2 will be unaffected by the index lowering operation ^ = 
dfivu" ') but its time component u° = dt/dr = 7 will differ in sign from the corresponding 
component uq = —7c 2 of the associated covector On the other hand the 3- momentum p a 
and energy E determine a 4- momentum covector n v with components 7r a = p a , /i = —E that 
are intrinsically covariant. The covariant nature of the momentum can be seen from the way 
it is introduced by the defining equation, 

5F- (48) 

in terms of the relevant position and velocity dependent Lagrangian function L, from which 
the corresponding equation of motion is obtained in the well known form 

dr dx» ■ 1 ' 

In the case of a free particle trajectory, and more generally for fluid flow trajectories in 
a model of the simple barotropic kind that is relevant in the zero temperature limit, the 
Lagrangian function will have the familiar standard form 

L = \^g^u^u u - , (50) 



in which (unlike what is needed for more complicated chemically inhomogeneous models [77, 
tf3j) it is the same scalar spacetime field fi that plays the role of mass in the first term and 
that provides the potential energy contribution in the second term. The momentum will thus 
be given by the simple proportionality relation 

\i v = [iu v , (51) 

so that one obtains the expressions E = 7/xc 2 , p a = fj,^/v a , in which the field fi is interpretable 
as the relevant effective mass. 

In the case of a free particle model, the effective mass /i will of course just be a constant, 
H = m. This means that if, as we have been supposing so far, the metric g^ u is that of flat 
Minkowski type, the resulting free particle trajectories will be obtainable trivially as straight 
lines. However the covarant form of the equations (^7j) to (^) means that they will still be 
valid for less trivial cases for which, instead of being flat, the metric g^ v is postulated to have 
a variable form in order to represent the effect of a gravitational field, such as that of a Kerr 



black hole (for which, as I showed in detail in a much earlier Les Houches school 78], the 
resulting non trivial geodesic equations still turn out to be exactly integrable). 
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In the case of the simple "barotropic" perfect fluid models with which we shall be concerned 
here, the effective mass field p will be generically non-uniform. In these models the equation 
of state giving the pressure P as a function of the mass density p can most conveniently be 
specified by first giving p in terms of the corresponding conserved number density n by an 
expression that will be decomposible in the form 

p = ran + 4j , (52) 
cr 

in which m is a fixed "rest mass" characterising the kind of particle (e.g. a Cooper type neutron 
pair) under consideration, while e represents an extra compression energy contribution. The 
pressure will then be obtainable using the well known formula 

P = (np- p)c 2 , (53) 

in which the effective dynamical mass is defined by 



dp 1 de 

a =-f=m+——. (54) 
an &■ an 



It is this same quantity p (sometimes known as the "specific enthalpy") that is to be taken 
as the effective mass function appearing in the specification ( |50|) of the relevant Lagrangian 
function (on what is formally identifiable as the tangent bundle of the spacetime manifold). 

When one is dealing not just with a single particle trajectory but a spacefilling fluid flow, 
it is possible and for many purposes desirable to convert the Lagrangian dynamical equation 



p9| ) from particle evolution equation to equivalent field evolution equations [77, [f^ ]. Since 
the momentum covector p u will be obtained as a field over spacetime, it will have a well 
defined gradient tensor V p p u that can be used to rewrite the right hand side of (^) in the 
form dp v /dr = u p V p p u . Since the value of the Lagrangian will also be obtained as a scalar 
spacetime field L, it will also have a well defined gradient which will evidently be given by an 
expression of the form 

dL dL „ 

We can thereby rewrite the Lagrangian dynamical equation ( [49| ) as a field equation of the 
form 

u p V p p v + p p V uU P = V V L . (55) 

An alternative approach is of course to start from the corresponding Hamiltonian function, 
as obtained in terms of the position and momentum variables (so that formally it should be 
considered as a function on the spacetime cotangent bundle) via the Legendre transformation 

H = p v u v - L . (56) 

In this approach the velocity vector is recovered using the formula 

dx^ dH 



dr dp u ' 

and the associated dynamical equation takes the form 

dp u dH 



(57) 



dr dx» ■ (58) 
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The consideration that we are concerned not just with a single trajectory but with a space- 
filling fluid means that, as in the case of the preceeding Lagrangian equations, so in a similar 
way this familiar Hamiltonian dynamical equation can also be converted to a field equation 
which takes the form 

2u p V [p » u] = -V U H, (59) 

with the usual convention that square brackets are used to indicate index antisymmetrisation. 
On contraction with u u the left hand side will evidently go out, leaving the condition 

u u V u H = , (60) 

expressing the conservation of the value of the Hamiltonian allong the flow lines. 

The actual form of the Hamiltonian function that is obtained from the simple barotropic 
kind of Lagrangian function (|50| ) with which we are concerned will evidently be given by 

H = ^ W + ^ , (61) 

in which it is again the same scalar spacetime field ^ that plays the role of mass in the first 
term and that provides potential energy contribution in the second term. 

In order to ensure the proper time normalisation for the parameter r the equations of 
motion (in whichever of the four equivalent forms (^9|), (|55|), (pgj), ( |59|) may be preferred) 
are to be solved subject to the constraint that the numerical value of the Hamiltonian should 
vanish, 

H = , (62) 

initially , and hence also by (^0|) at all other times. This is evidently equivalent to imposing 
the standard normalisation condition 

u% = -c 2 , (63) 
on the velocity four vector. In more general "non-barotropic" systems, such as are needed 



for some purposes, the Hamiltonian may be constrained in a non uniform manner [77, |73| so 
that the term on the right of ( |59| ) will be non zero, but in the simpler systems that suffice 
for our present purpose the restraint ( |62|) ensures that this final term will drop out, leaving a 
Hamiltonian equation that takes the very elegant and convenient form 

u"w vp = . (64) 

in terms of the relativistic vorticity tensor that is defined as the antisymmetrised ( "exterior" ) 
derivative of the momentum covector, i.e. 

w vp = 2V^/x p] . (65) 

It is an evident consequence (and, as discussed in greated detail in the above cited Les 
Houches notes [55], would still be true even if ( |62| ) were not satisfied) that if is zero 



initially it will remain zero throughout the flow, which in this case will be describable as 
"irrotational" . 

In cases for which the vorticity is non-zero, the "barotropic" dynamical equation ( |64|) is 
interpretable as requiring the flow vector vP" to be a zero eigenvalue eigenvector of the vorticity 
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tensor w„ vi which is evidently possible only if its determinant vanishes, a requirement that is 
expressible as the degeneracy condition 

«V[i/*V] = • ( 66 ) 

Since the the possibility of it having matrix rank 4 is thus excluded, it follows that unless it 
actually vanishes the vorticity tensor must have rank 2 (since an antisymmetric matrix can 
never have odd integer rank). This means that the flow vector is just a particular case 
within a whole 2-dimensional tangent subspace of zero eigenvalue vorticity eigenvectors, which 
(by a well known theorem of differential form theory) will mesh together to form well defined 
vorticity 2-surfaces as a consequence of the Poincare closure property, 

V[n w M = > ( 67 ) 

that follows from the definition (|65[). 

Although it has long been well known to specialists p^j , the simple form (34) of what is 



interpretable just as the relativistic version of the classical Euler equation is still not as widely 
familiar as it ought to be, perhaps because its Hamiltonian interpretation was not recognised 



until relatively recently 1 77], . It does not constitute by itself the complete set of dynamical 
equations of motion for the perfect fluid, but must be supplemented by a particle conservation 
equation of the usual form for the particle number current 

n v = nu v , (68) 

which must of course satisfy the condition 

V u n" = . (69) 

A much more widely known, but for computational purposes (particularly in curved space- 
time) less useful form of the perfect fluid dynamical equations is to express them in terms of 
the stress momentum energy density tensor, which is given in terms of the mass density p and 
the pressure P by 

= (p+ ^)u»u u + Ptf v , (70) 
cr 

and which must satisfy a so called conservation law of the standard form 

y u T^ v = . (71) 

Although it is conveniently succinct, a disadvantage of this traditional formulation is that it 
is directly interpretable as a law of conservation of momentum and energy in the strict sense 
only in the case of a flat (Minkowski type) spacetime, but not in a curved background such 
as that of a neutron star. The possibility in the barotropic case (i.e. when P is a function 
only of p) of decomposing the combined set of dynamical equations (|70| ) as the combination 
of the convergence condition (|69|) (obtained by contracting (^TJ) with Up) and the relativistic 
Euler equation (|64]), which can be written out more explicitly as 



n»V [u p p] =0, (72) 

has the advantage that these are interpretable as genuine conservation laws - for particle 
number flux and vorticity respectively - even in an arbitrarily curved spacetime background. 
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3 Single constituent superfluid models 



The simplest superfluid models, namely those pertaining to the zero temperature limit, are 
just ordinary perfect fluid models subject to the restraint of irrotationality, with a momentum 
covector given as the gradient 

Hu = VuS, (73) 

of a scalar field S. This scalar field is to interpreted as being proportional to the angle of 
the mesoscopic phase factor, e 1 ^ say, of an underlying scalar bosonic condensate, in which the 
phase angle eft is given according to the usual correspondence principle by 

</> = S/h. (74) 

In the most familiar application the bosons are Helium-4 atoms, while between the ions of a 
neutron star crust below the neutron drip transition they will be Cooper type neutron pairs. 
(However a less simple description is not sufficient for the spin 1 neutron pairs below the base 



of the crust, nor in the even more complicated, though experimentally accessible, case [67] of 
Helium-3, for which a relativistic description is still not available). 

In a multiconnected configuration of a classical irrotational fluid the Jacobi action field 
S obtained from ( fTq) might have an arbitrary periodicity, but in a superfluid there will be a 
U(l) periodicity quantisation requirement that the periodicity of the phase angle cj> should be 
a multiple of 2tt, and thus that the periodicity of the Jacobi action S should be a multiple of 
2nh. The simplest configuration for any such superfluid is a uniform stationary state in a flat 
Minkowski background, for which the phase will have the standard plane wave form 

S/h = k a x a - cot , (75) 

from which one obtains the correspondence \x v <-> {—Tito, hk a }, which means that the effective 
energy per particle will be given by E = 7/ic 2 = hu> and that the 3-momentum will be given 
by p a = n-yva = hk a . 

It is to be remarked that for ordinary timelike superfluid particle trajectories the corre- 
sponding phase speed u/k of the wave characterised by ( |75|) will always be superluminal - a 
fact of which people working with laboratory Helium-4 tend to be blissfully unaware, and can 
usually safely ignore, since what matters for most practical purposes is not the phase speed 
but the group velocity of perturbation wave packets. 

In the irrotational case characterised by (^) the Euler equation (|64|) is satisfied automat- 
ically, so the only dynamical equation that remains is (^). When the phase scalar is subject 
to a small perturbation, 8<f> = (p say, it can be seen that the corresponding perturbation of 
the conservation law ([39]) provides a wave equation of the form 

D<P = 0, (76) 

in which □ is a modified Dalembertian type operator that is constructed from an appropri- 
ately modified space-time metric tensor g^ v in the same way that the ordinary Dalembertian 
operator □ = V^V^ is constructed from the ordinary spacetime metric tensor g^ v '. The ap- 
propriately modified spacetime metric, namely the relativistic version of what is known in the 
context of Newtonian fluid [|79|, ISO] and superfluid |8l]] mechanics as the Unruh metric, can be 



read out in terms of the light speed c and the (first) sound speed c : given by fl46|) as 



= eff- c^W), (77) 
n 
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where 7^ is the spacially projected (positive indefinite) part of the ordinary space time metric, 
as defined by 



The quantum excitations of the linearised perturbation field cp governed by (|7q ) are what 
are known as phonons. For such excitations the phase speed and the group velocity are the 
same, both being given with respect to the unperturbed background by the ordinary ("first") 
soundspeed, c t , as given by (|46|), which will of course be subluminal. Phonons do nevertheless 
have a tachyonic aspect of their own, because the fact that their phase speed is subluminal 
automatically implies that they have a 4-momentum covector that is spacelike, in contrast 
with that of a ordinary fluid or superfluid particle which is timelike. This means that whereas 
the effective energy E of an ordinary fluid or superfluid particle is always positive, the effective 
energy E of a phonon may be positive or negative, depending on whether the frame of reference 
with respect to which it is measured is moving subsonically or supersonically. The well known 
implication is that if the superfluid is in contact with a supersonically moving boundary there 
will inevitably be an instability giving rise to dissipative phonon creation. 

Given a dynamical system, one of the first things any physicist is inclined to ask is whether 
it is derivable from a Lagrangian type variation principle. We have already seen in the previous 
sections that (^) by itself is obtainable from Lagrangian equations of motion for the individual 
trajectories, which are of course obtainable from a one dimensional action integral of the form 
/ Ldr with L as given by ([)(]). The question to be adressed now is how to obtain the complete 
set of dynamical equations (|7l[), including ( p4| ) as well as (|69|), from an action integral over 
the 4-dimensional background manifold of the form 

CdS^, dS& = Mll d * Xj (79) 

c 

for some suitable scalar Lagrangian functional C. 

There are several available procedures for doing this for a generic perfect fluid with rota- 
tion, involving radically different choices of the independent variables to be varied: although 
they are all ultimately equivalent "on shell" the "off shell" bundles over which the variations 
are taken differ not only in structure but even in dimension. These methods notably include 
the worldline variation procedure (the most economical from a dimensional point of view) de- 



veloped by Taub[83], and the Clebsch type variation proceedure developed by Schutzp4[, as 
well as the more recently developed Kalb-Ramond type method []85|, [55| that has been specif- 
ically designed for generalisation (54|] to allow for the anisotropy arising from the averaged 
effect of vortex tension in the treatment of superfluidity at a macroscopic level. None of these 
various methods is sufficiently simple to have become widely popular. 

The problem is much easier to deal with if, to start off with, one restricts oneself to the 
purely irrotational case characterised by (|7^), which is all that is needed for the description 
of zero temperature superfluidity at a mesoscopic level. In this case a very simple and well 
known procedure is available. In this procedure, the independent variable is taken to be just 
the Jacobi action S, or equivalently in a superfluid context, the phase (p as given by ([74]), and 
the action is simply taken to be the pressure P expressed as a function of the effective mass ji, 
with the latter constructed as proportional to the amplitude of the 4-momentum, according 
to the prescription 

M V = -/j^n" , (80) 
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with the 4- momentum itself given by the relation (73) that applies in the irrotational case, 
i.e. 

H v = hV v <f>. (81) 

Thus setting 

C = P , (82) 
and using the standard pressure variation formula 

SP = c 2 n5p , (83) 

one sees that the required variation of the Lagrangian will be given by 

5C = -n v 6fi v = -hn v V v {54>) . (84) 

Demanding that the action integral ( |79|) be invariant with respect to infinitesimal variations 
of (p = Sep then evidently leads to the required conservation law fl69|). 

It is to be noted that this variational principle can be reformulated in terms of an in- 
dependently variable auxiliary field amplitude <3? and an appropriately constructed potential 
function V^j^} as a function of which the action takes the desirably fashionable form 

-C = -Y* 2 (V^)V>-F{*}, (85) 

which is interpretable as the classical limit of a generalised Landau Ginzburg type model. In 
this formulation, as discussed in greater detail elsewhere [85, 55], the auxiliary amplitude is 
to be identified as being given by the formula 



y/p + Pjc 1 V 



(-) , (86) 



while the prescription for the corresponding potential energy density function is that it should 
be given by 

V = ^- (87) 

Having evaluated V as a function of $ one can recover the effective mass //, number density 
n, mass density p and pressure P of the fluid using the formulae 

2 1 dV o 

and 



p=i$V + -^, P =\<S> 2 i?c 2 - V . 



which are derivable from (|53| ) and (M). It is to be remarked that the covariant inverse of the 
generalised Unruh tensor ( |77| ) is expressible in this notation as 

<T M l = ^(Cj'V* - c lC - 2 u^u u ) . (90) 



A particularly noteworthy example is the conformally invariant special case [55] char- 
acterised by a potential function that is homogeniously quartic, V oc $ 4 , which is what is 
obtained for a radiation gas type equation of state of the familiar form P = pc 2 /3, and for 
which the (first) sound speed is given by c 2 = c 2 /3. 
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4 Landau type 2-constituent superfluid models. 



As an intermediate step between the very simple single constituent superfluid models described 
in the previous section and the more elaborate models needed in the context of neutron stars, 
the purpose of this section is to describe the relativistic version of the category of non dissi- 
pative 2-constituent superfluid that was originally developed by Landau for the description of 
ordinary superfluid Helium-4 at non-zero temperature. As well as the relevant conserved par- 
ticle number current (representing the flux of Helium atoms in that particular application) 
such a model involves another independently conserved current vector, say, representing 



the flux of entropy. In the single constituent case characterised by the variation rule (84) we 
saw how the current vector n u was associated with a dynamically conjugate covector fi u that 
is interpretable as representing the effective mean 4- momentum per particle. In a similar way 
in a 2-constituent model the second current vector s u will be analogously associated with its 
own dynamically conjugate 4-momentum covector B^. 

The earliest presentations of the generic category of non-dissipative 2-constituent super- 
fluid were on one hand a generalisation [[74|] of the relativistic Clebsch formulation |84| based 
on the variation of a generalised pressure function ^ depending on the 4-momentum covectors 
fjL u and according to the partial differentiation rule 

d$> = - n v dfi u - s u d@ u , (91) 

and on the other hand a generalisation of the world line variational formulation due 



to Taub [83] based on the variation of a master function A depending on the currents and 



s v according to the partial differentiation rule 

dA = jx v dn v + S v ds v . (92) 

Although they were originally developed independently these alternative formulations were 
subsequently shown to be equivalent to each other and to an intermediate crossbred ver- 
sion [|rj| based on a Lagrangian density 

C = ¥ + s u Q u = A - n^fiu , (93) 

depending on the particle 4-momentum covector \i v and the entropy current s u according to 
the partial differentiation rule 

dC = Q u ds v - n v dfi u . (94) 

All of these variational formulations are subject to the complication that the allowable 
field variations are not free but must be suitably constrained to avoid giving overdetermined 
field equations. Although it violates the symmetry between the two kinds of conserved current 
n v and s u that are involved, the crossbred formulation characterised by (^) is the one that 
allows the simplest specification of the constraints required to get the appropriate dynamical 
equations for the superfluid case. In this formulation |75| the constraint on the particle 4- 
momentum covector is simply that it should have the same phase gradient form ([T3]) as in the 
zero temperature limit in which the entropy constituent is absent, namely 

\t v = hV v <f> . (95) 

The corresponding constraint on the current vector s v of the "normal" constituent is the 
not quite so simple Taub type requirement that its variation should be determined by the 
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displacement of the flow lines generated by an arbitrary vector field C, v say, which means [73] 
that it must have the form 

ds v = ( p V p s u - s v V p C + s u V p ( u , (96) 

whose derivation is obtainable by a procedure that will be explained more explicitly in the 
next Section. Demanding invariance of the volume integral of C with respect to infinitesimal 
local variations of the phase variable (ft then gives the usual particle conservation law in the 
same form (|69"| ) as for the single constituent limit, while demanding invariance for an arbitrary 
local displacement field £ gives not only the analogous entropy conservation law 

V u s u = , (97) 

but also the dynamical equation, 

s u V [u Q p] =0, (98) 

that governs the evolution of the thermal 4-momentum covector in a manner analogous to 
that whereby the relativistic Euler equation (|7^) governs the evolution of the momentum 
covector in an ordinary perfect fluid. These dynamical equations entail (but unlike the single 
constituent case are not entirely contained in) an energy momentum pseudo-conservation law 
of the usual form (fn|) for a stress-momentum-energy density tensor that can be written in 
the form 

Tp = ny p + s»Q p + ^g v p , (99) 

which will in fact (although it is not obvious in this particular expression) be automatically 
symmetric, = 0. 

The category of models characterised by the preceeding specifications for various conceiv- 
able forms of the equation of state specifying £ as a scalar function of [i v and s u is very large. 
The use of what is interpretable p7[ as a special subcategory therein, on the basis of a partic- 



ular kind of separation ansatz, was proposed in early work of Israel [38] and Dixon [89] and 
has been advocated more recently by Olsen |Xj]. Unfortunately however, the simplification 
provided by the Israel Dixon ansatz (effectively the relativistic generalisation of the obsolete 
Tisza-London theory that was superceded by that of Landau) is incompatible with the kind of 
equation of state that is needed for even a minimally realistic treatment of a real superfluid. 

A satisfactory treatment of what goes on at temperatures high enough for non-linear "ro- 
ton" type excitations to be important is not yet available, but in the low temperature "cool" 
regime, in which only linear "phonon" type excitations are important, it is not difficult to 
provide a straightforward analytic derivation of the kind of equation of state that is appropri- 
ate. Following the lines developed in a non - relativistic context by Landau himself |)l[] the 
relativistic version of the approriate "cool" equation of state has recently been derived [Jf™ by 
considering perturbations of the single constituent model - with equation of state specified as 
a pressure function, P{/j-} ~ that describes the relevant zero temperature limit. The result is 
obtained in an analytically explicit form that (despite the fact that it is not of the separable 
Israel Dixon kind) can be given a very simple expression in terms of what we referred to as 
the "sonic" metric, which is specifiable by the conformal relation 



in terms the Unruh phonon metric (77) that is associated with the relevant zero temperature 



limit state as specified by the relevant momentum covector \i v , which by (pQ) determines the 
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relevant value of the scalar p and hence (via the zero temperature equation of state, using the 
formalism of Section |2.2| ) also of the relevant phonon speed c t and field amplitude <£. While 
the Unruh metric is more convenient for many purposes, the advantage of the conformal 
modification we have used, namely 

gP° = gP° + (c~ 2 - c- 2 )u p u a , (101) 

is that its spacially projected part agrees with that of the ordinary space metric, from which 
it differs only in the measurement of time. 

The result that is obtained [76] is given by a Lagrangian of the form 

C = P-3ip (102) 

in which the deviation from the zero pressure limit value P{p} is given as a function not just 
of the particle 4-momentum covector p u but also of the entropy flux s v (postulated to be 
sufficiently weak to be constituted only of phonons) by the formula 

^ = \c~^\g-y s ^i\ (103) 

where fi is identifiable to a very good approximation with the usual Dirac-Planck constant h, 
its exact value being given by 

*=£(fr =">■»»■ < im > 

This is equivalent to taking the generalised pressure function to be 

$ = P + <0, (105) 

with 

V> = f (J^V^e,) 2 , (106) 

in which the effective thermal 4-momentum per unit of entropy is given (according to the 
partial differentiation formula (|94|) by 



dti 

with 

c 2 

Q- P l = g P o + (i - -^r) i^vrVpA** • ( 108 ) 

An concrete illustration, allowing the explicit evaluation of the relevant quantities, is 
provided by the polytropic case, as characterised by a (single consitituent) equation of state 
giving the mass density p as a function of the number density n in terms of a fixed ("rest") 
mas per particle m, a scale constant k and a fixed dimensionless index 7 in the form 

p = ran + ktl 1 44> p = m + Kjn' y ~ 1 , (109) 

which corresponds to taking the pressure to be given by 

P = kc 2 ( 7 - 1W = *c 2 ( 7 - l)f V—!^) 7 ^ , (HO) 

V K7 / 

while the corresponding sound speed will be given (independently of k) by 

Cl 2 = ( 7 -l)(l-^)c 2 . (Ill) 
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5 Non- conservative model with transfusion and vortex drag 



Although the Landau type of model described in the previous section has been found to be 
very effective for the description of liquid Helium-4 under laboratory conditions, it is not of 
much use for direct application in neutron star matter because the thermal effects it allows 
for will in general be less important than other complications whose treatment will require 
the use of more elaborate models whose relativistic versions are still at a relatively early stage 
of development and will not be presented here. The most important of these complications, 
whose treatment in a relativistic framework has been the subject of preliminary work that is 
discussed elsewhere, are due to the effect of the protons that will be present, either in ionic 
nuclei that are responsible for the elastic solid behavior |T3| of the crust, or as a dissolved 



superfluid [[7C], 71 1 at deeper levels. Another complication that is relevant for the macroscopic 
treatment of a neutron star is the necessity of averaging over an Abrikosov type lattice of 
quantised vortices (that must be roughly aligned with the rotation axis of the star) whose 
effective tension entails deviations |54] , [55| from perfect fluid isotropy. 

Like the thermal effect discussed in the preceding section, these various complications 
can all be provisionally set aside as perturbations to be incorporated at a later stage in a 
systematic approach whose first stage requires the use only of a relatively crude description 
in which, except for the superfluid neutrons with baryon number current vector n„ all the 
other constituents, meaning mainly protons and electrons, move together with the entropy 
as a single "normal" constituent with baryon number current n". Whereas the anisotropy 



arising from vortex tension [54, [55[ is relatively unimportant, a major role in the long term 
evolution of the star is likely to be played by the static pinning or dynamical drag forces 
exerted on the vortices by the composite "normal" background constituent. Another effect 
that is of importance in the long run is that of "transfusion" whereby - due to the subduction 
resulting from the drift mechanism whose effect is roug hly described by © or @ - the 
superfluid neutron contribution to the baryon current may undergoes transformation (via 
weak beta decay type processes) to the "normal" (essentially protonic) consituent, and vice 
versa, so that only the total baryon current 

n b = < + < (112) 

remains locally conserved throughout: 

V„n£ = 0. (113) 

The kind of (non-conservative) 2-constituent model needed for this purpose is obtainable 
as a generalisation of the kind of (conservative) 2-constituent superfluid model discussed in 
the preceding section, starting from the formulation in terms of a master function A in which 
the currents (not momenta) are taken as the independent (but not entirely free) variables. 

In a transfusive model of the type set up here, the "normal" constituent is not entirely 
dependent on (though it does include) entropy, so that it is present even at zero temperature: 
the primary role of this non - superfluid constituent is to represent the fraction of the baryonic 
material of the neutron star that is not included in the neutron superfluid, as well as the 
degenerate electron gas that will be present to neutralise the charge density resulting from 
the fact that some of these baryons will have the form of protons rather than neutrons. In 
the solid "crust" layers of a neutron star the protons will be concentrated together with a 
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certain fraction of the neutrons in discrete nuclear type ions, which at the relatively moderate 
temperatures that are expected to apply will form a solid lattice. In the upper crust the 
"normal" constituent consisting of the ionic lattice and the degenerate electrons will include 
everything, but in the lower crust (at densities above about 10 11 gm/cm 3 ) the crust will 
be interpenetrated by an independently moving neutron superfluid. What we refer to as 
"transfusion" occurs when compression takes place so that the ionic constituent undergoes a 
fusion process whereby neutrons are released in the form of newly created superfluid matter, 
or conversely, when relaxation of the pressure allows excess neutrons to be reabsorbed into 
the ions. 

A more elaborate treatment would specifically allow for the expectation that the pro- 
tons would form an independently conducting superfluid of their own at very high densities, 
whereas they will combine with some of the neutrons at intermediate densities, and with all 
of the neutrons at low densities, to form discrete ions which will tend to crystalise to form 
a possibly anisotropic lattice. What matters for our present purpose is that regardless of its 
detailed constitution, all this "normal" matter will in effect be strongly self coupled (T^j by 
short range electromagnetic interactions so that its movement will be describable to a very 
good approximation as that of a single fluid with a well defined 4-velocity, u c ^ say, the only 
independent motion being that of the (electromagnetically neutral) neutron superfluid with 
velocity u n ^ say. The latter will specify the direction of the part of the baryon current, 

= n n u^ , (114) 

carried by the neutron superfluid, while the "normal" matter velocity specifies the direction 
of the remaining collectively comoving part, 

n/ = n c u/, (115) 

of the baryon current. 

At densities below the "neutron drip" transition at about 10 11 gm/cm 3 , the "normal" 
collectively comoving constituent n c ^ will of course be identifiable with the total, . The 
reason why the remaining free neutron part n n M - which will always be present at higher 
densities - is presumed to be in a state of superfluidity is that the relevant condensation 
temperature, below which the neutrons form bosonic condensate of Cooper type pairs is 
estimated [92] to be at least of the order of 10 9 K, while it is expected that a newly formed 
neutron star will drop substantially below this temperature within a few months [p3|]. At such 
comparatively low temperatures the corresponding entropy current s M say will not play a very 
important dynamical role, but for the sake of exact internal consistency it will be allowed 
for in the model set up here, in which it will be taken for granted that it forms part of the 
"normal" collectively comoving constituent so that it will have the form 

= s< . (116) 

Under conditions of sufficiently slow convection, the transfer needs not involve significant 
dissipation, so the process should be describable by a Lagrangian scalar, A say, that will 
depend just on the currents introduced above, of which the independent components are given 
just by the vectors n£ and and the scalar s. As a first approximation (whose accuracy in 
the various relevant density regimes is a subject that needs much further investigation) one 
might suppose that the Lagrangian separates in the form A = —p c c 2 — p n c 2 in which p c is a 
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mass density depending only on s and n c , while p n is an another energy mass depending only 
on n n , but we shall not invoke such a postulate here, i.e. we allow for the likelihood that the 
properties of "normal" constituent will be affected by the presence of the superfluid constituent 



and vice versa, which means that there will be an entrainment effect p4|, [lq , 19, [95|], whereby 
for example the velocity of the superfluid neutron current will no longer be parallel to the 
corresponding momentum. (As an alternative to the more suitable term "entrainment" this 
mechanism is sometimes referred to in the litterature as "drag" , which is misleading because 
entrainment is a purely conservative, entirely non-dissipative effect, whereas the usual kinds 
of drag in physics, and in particular the kind of drag to be discussed below, are essentially 
dissipative processes.) 

If we adopted the (gas type) description embodied in the separation ansatz we would 
have two separate variation laws which in a fixed background would take the form c 2 5p c = 
<B5s + c 2 p c 5n c and c 2 5p n = c 2 p n Sn n , in which would be interpretable as the temperature, 
/i° would be interpretable as the effective mass per baryon in the "normal" part, and p, n would 
be efective mass per neutron in the superfluid part (which would be equal to its analogue in 
the "normal" part, i.e. p n = p c , in the particular case of a state of static thermodynamic 
equilibrium.) 

In the less specialised (liquid type) description to be used here, there will just be a single 
"conglomerated" variation law, whose most general form, including allowance for a conceivable 
variation of the background metric, will be expressible (correcting one of the copying errors 
in the originally published version j5(|) as 

5A = -Q5s + p\bn^ + /A<5< + i(cr 2 0s« + n»p cv + n n ^ nu )S 9flu , (117) 

where is to be interpreted as the temperature and where and p ^ are to be interpreted 
as the 4-momentum per baryon of the neutron superfluid and the "normal" constituent re- 
spectively. 

To obtain suitable fluid type dynamical equations from a Lagrangian expressed as above 
just in terms of the relevant currents, the variation of the latter must be appropriately con- 



strained in the manner[73] that was originally introduced for the case of a simple perfect fluid 
by Taub[83|. The standard Taub procedure can be characterised as the requirement that 
the variation of the relevant current three form, which for the "normal" constitituent in the 
present application will be given in terms of the antisymetric space-time measure tensor e^pa 
by ' ' 

N/jjsp = E^vpcjUl , (118) 

should be given by Lie transportation with respect to an associated, freely chosen, displace- 
ment vector field (J* say. This ansatz gives the well known result 

SN^p = C A V A AW + 3iV A ^V p] C A • (119) 

Although a variation Sg^ of the metric has no effect on the fundamental current three form, 
Nfjivp, it will contribute to the variation of the corresponding vector, 

< = ^P°N up(T , (120) 

for which one obtains 

Sri? = - <V„C + <(V„<r - h up 5g up ) . (121) 
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(Application of an analogous procedure to the entropy current provides the variation rule 
that was used in Section ||).) In terms of the orthogonally projected metric, 

7c ^ = ^ + c- 2 <<, (122) 
the corresponding variation of the unit flow vector will be given by 

6uS = 1c %(CV uU p - U( yv u C p ) + ^c- 2 u^ c u p c 5 gi/p , (123) 
and the corresponding variation in the current amplitude n c will be 

5n c = V.KC) + n c (c~ VW M C„ - \l^ v • (124) 

Since the entropy flux is to be considered as comoving with the "normal" constituent, it is 
subject to a variation given by the same displacement vector which thus gives 

5s = V„(aO + s(c~V W^G - \ir^u) ■ (125) 

On the other hand for the superfluid constituent there will be an independent displacement 
vector field (J* say, in terms of which the analogously constructed variation will be 

k = rv,< - Kv„e + <(v,r - l<f p s gvp ) . (126) 

The effect of this variation process on the Lagrangian density HgH^A itself can be seen 
to be expressible in the standard form 

iisir 1/2 *([i0ii 1/2 A) = cfu + er* + + , (127) 

in which f c u will be interpretable as the force density acting on the "normal" constituent, f n u 
will be interpretable as the force density acting on the superfluid constituent, and T^" will 
be interpretable as the stress momentum energy density of the two constituent as a whole. 
By considering the trivial case in which there is no actual physical alteration of the system, 
but in which the apparent changes are merely due to the displacement of the reference system 
generated by a vector field £ v = ( u , in which case the apparent variation of the metric will be 
given by 5g uu = 2V^(^, it can be seen from ( 127 ) that the separate forces must automatically 
satisfy an identity of the form 

f v +f\ = r\, (128) 

where f e \ is the conglomerated external force density that is defined by 

! c \ = V M T^ . (129) 

The residual current TV 1 in the divergence will be of no importance for our present purpose 
(by Green's theorem it just gives a surface contribution that will vanish by the variational 
boundary conditions) but it is to be noted for the record that it will have the form 

TV = 2(K c ^(c- 2 9su cl/ + n c /4) + 2£ ^n n V, . (130) 

The conglomerated stress momentum energy density tensor can easily be read out as 

T" u = %t + c~ 2 ®suSu cv + ra c V*/ + n^^ v , (131) 
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where 

* = A + se - n c v fi c v - n n V n , • (132) 

(Although this expression is not manifestly symmetric, the asymmetric contributions will 
automatically cancel due to the identity What matters most for our 

present purpose is the form of the respective force densities: the force law (i.e. the relevant 
relativistic generalisation of Newton's "second" law of motion) for the "normal" constituent 
is found to take the form 

f% = 2s^{c~ 2 ®u cu] ) + 2n c ^fi c v] + c- 2 Qu cu V^ + fi c u V^ , (133) 

while the force law for the superfluid component is found to take the simpler form 

fl = ft + fT, (134) 

in which the first term is a "chemical" contribution, representing the effect of any neutron 
superfluid particle creation or destruction, which is given by 

A=^VX. (135) 



The last term in (134) is a "mechanical" contribution, allowing for drag or pinning forces 
exerted on the vortices by the crust and balanced by the Magnus effect, according to the 
formula 

f me , = n^w%, (136) 

using the notation 

w% = 2V [M M*] (137) 

for the vorticity 2-form of the superfluid neutrons. It is to be noted that this is not the 
mesoscopic (intervortex) superfluid vorticity, which simply vanishes, but the average vorticity 
on a macroscopic scale that is large compared with the spacing (typically a very small fraction 
of a cm.) between the superfluid vortices. For a very accurate treatment it would be necessary 
to take account of the macroscopic anisotropy resulting from the effective tension of these 
vortices, as has already been done [54, [5^] for the case a single constituent, but for the 



discussion of global evolution on timescales long compared with the stellar oscillation periods 
(a small fraction of a second) such an effect seems unlikely to be important. 

Although the complete expression (|133| ) is not so simple, it is to be observed that the 
time component in the "normal" rest frame (representing the rate of working on the "normal" 
constituent) as obtained by contraction with the relevant unit vector u c u has the comparitively 
simple form 

u c v f c v = u c v n c u V^nS - 0V M s" . (138) 

If we were to impose the variation principle to the effect that the system should be invariant 
with respect to arbitrary worldline displacements (as specified by the independent fields Q u 
and it would follow that each of the forces /° and /° would have to vanish. However it is 
evident from the identity ( |128| ) that we cannot adopt such a restrictive postulate in a model 
designed to treat the effect of pulsar slowdown due to a torque attributable to coupling to an 
external electromagnetic field that is removing angular momentum by radiation to infinity. 
As well as the intrinsically non-conservative magnetic torque contribution to it is also 
important Eft to include a contribution to allow for the effect of the elastic solidity in the 
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crust, which is not incorporated into the simple fluid type model included here (and which 
would require the use of a much more elaborate model [73] for its detailed evaluation). 



Although our ultimate purpose is to allow for a non vanishing external torque force, 
whatever force law we assume must be such that if the external force were somehow 
switched off so as to leave an effectively isolated system, the second law of thermodynics (no 
decrease of entropy in an isolated system) would be respected, i.e. we must have /® x = => 
V u d u > 0. It can be seen that this is equivalent to the requirement of positivity of the right 
hand side of the identity 

0V ^ + U( y ru = U( y^ u _ ^jv^nn" + u c v r\ , (139) 



that is obtained from ( |138| ), taking account of the total baryon conservation law (|112| ). Since 
they involve very different physical processes, one comes to the conclusion that each of the 
two terms on the right of (|13S| ) must satisfy its own separate positivity condition. 

The positivity requirement for the first of these terms is presumably to be attributed to a 
crust particle creation law of the form 

V M < = H<(/x\ - n%) (140) 

for some positive coefficient S. Such a law is an obviously natural generalisation of the kind 
of creation rate formula that is familiar in chemical physics. In the present context what 
is involved is conversion of protons to neutrons by weak interactions, and the situation is 
complicated by the consideration that as far as the large scale mechanics of the neutron star 
is concerned, the effective rate may depend not just on microscopic processes, but also, when 
subduction is involved, on the rather messy process whereby the crust is broken up before it 
ultimately dissolves. 

To complete the specification of the system, all that remains is to find the appropriate 
ansatz for the mechanical force f me v . This problem is more delicate than that of the (effec- 
tively scalar) chemical case, since as well as the "second law" requirement u"f me u > 0, the 
answer must respect the nature of the macroscopic vorticity 2-form w n ^ v which although non 
vanishing (unlike the mesoscopic vorticity between vortices) cannot be arbitrary (as in an 
ordinary viscous fluid): to be consistent with the underlying superfluid nature of the neutron 
constituent, it must satisfy an algebraic degeneracy condition of the form ( |66|) in order to be 
compatible with the existence of a well defined congruence of orthogonal 2-surfaces generated 
by (non vanishing) tangent vectors, v v say, such that w ri )ll ,v u = 0. It can be seen from the 
form of the defining relation ( |136| ) that the obvious way to obtain this degeneracy property 
is to take the force law to have the form / me t/ = w„ a v a for some suitably chosen vector v v 
which, to satisfy the "second law" requirement must satisfy u^w I l IJ v a > 0. The required 
ansatz can thus be taken to be given by = aw^u" for some positive coefficient a. This 
result is conveniently expressible in terms of the rank-2 tensor lt l u of orthogonal projection 
with respect to the vortex 2-surface, which is given by 

±t= 2(w^ W l a )- 1 w^w\ u . (141) 

We end up with an expression taking the form 

r% = th ± vo < , (142) 
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for a positive resistive drag coefficient r/ r (given in terms of the previous coefficient a by 2r/ r 
= aw n P a w n pa ). 

The generic class of dissipative models characterised by finite values of H and r) T has 
four different kinds of non dissipative limit. In the low reactivity limit H — » we have the 
non-transfusive limit characterised by the separate superfluid particle conservation law 

V„< = , (143) 

whereas in the opposite high reactivity limit S — > oo we have the chemical equilibrium limit 
characterised by 

<(^"^) = 0. (144) 

which is what would be expected in cases for which the (continental drift like) crust circulation 
responsible for the transfusion is characterised by timecales that are very long compared with 
those p8| of the relevant weak (direct or inverse beta decay) interactions. For each of these 
conceivably relevant possibilities, we have the drag free limit tj t — > characterised by the 
condition of vanishing Magnus force, 

<w% = 0, (145) 

or at the opposite extreme the perfect vortex pinning limit, r) x — > oo characterised by the 
condition that the vortex worldsheets should be at rest with respect to the "normal" (crust) 
background, 

«, = 0. (146) 

The actual evaluation of the drag coefficient rj T , and the question of whether one or other 
of these simple extreme limits is realistic depends on delicate technical issues |96|, |9^] whose 
definitive resolution is not entirely clear. A scenario of the type envisaged by Ruderman |39|1 , 
as represented by Figure 1, and described by (^) (which seems to be appropriate for Vela) 
is what would be obtained in the case characterised by ( |146j ), whereas the more recently 
proposed alternative scenario |4(| represented by Figure 2, and described by ([lC|), is what 



would be obtained in the case characterised by ( 145 ). 

It is of course to be expected that such extreme scenarios will turn out in practise to 
be oversimplifications of a more complicated reality, whose description is likely to require 
modelisation with not just two |H] but many independently rotating components, to allow 
for the variation of the chemical and mechanical coefficients E and rj r over a wide range of 
finite values function of depth in the star. 

Acknowledgements The author wishes to thank Silvano Bonazzola, David Langlois, 
Eric Gourghoulon, Pawel Haensel, Isaac Khalatnikov, Reinhardt Prix, and David Sedrakian 
for conversations and collaboration. 



36 



References 

B.K. Harrison, K.S. Thorne, M. Wakano, J. A. Wheeler, Gravitation theory and gravita- 
tional collapse (University of Chicago Press, Chicago, 1965). 

A. B. Migdal, "Superfluidity of the moments of inertia of nuclei", Nucl. Phys. 13, 655-674 
(1959). 

N. K. Glendenning, Compact stars (Springer, Heidelberg, 1997). 

F. Weber, "Pulsars as astrophysical laboratories for nuclear and particle physica" (I.O.P. 
Publishing, Bristol, 1999). 

D. Langlois, "Superfluidity in relativistic neutron stars". [ astro-ph /0008 1 6 1 ] 

L.D. Landau, E.M. Lifshitz, Course of Theorteical Physics, Vols. 1-9 (Pergamon, Oxford, 
1959). 

G. Baym, D. Pines, "Neutron starquakes and pulsar speedup", Annals of Physics 66, 
816-835 (1971). 

G. Baym, C. Pethick, D. Pines, M. Ruderman, "Spin up in neutron stars: the future of 
the Vela pulsar", Nature 224, 872-874 (1969). 

M. Hoffberg, A.E. Glassgold, RW. Richardson, M. Ruderman, "Anisotropic superfluidity 
in neutron star matter", Phys. Rev. Lett. 24, 775-777 (1970). 

J. A. Sauls, D.L. Stein, JW. Serene, "Magnetic vortices in a 3 P2 neutron superfmid", 
Phys. Rev. D25, 967-975 (1982). 

N.K. Glendenning, "First order phase transitions with more than one conserved charge: 
consequences for neutron stars", Phys. Rev. D46, 1274-1287 (1992). 

N.K. Glendenning, S. Pei, '"Crystalline structure of the mixed confined - deconfined 
phase in neutron stars", Phys. Rev. C52, 2250-2253 (1995). 

M. Alford, K. Rajagopal, F. Wilczek, "Color-flavor locking and chiral symmetry breaking 
in high density QCD", Nucl. Phys. B537, 443-458 (1999). [|hep-th/9804403]| 

D. Blaschke, D.M. Sedrakian, K.M. Shahabasyan, "Diquark condensates and the mag- 
netic field of pulsars", Astron. Astrophys. 350, L47 (1999). [|astro-ph/9904395|1 

M. Alford, J. Berges, K. Ragjagopal, "Magnetic fields within colour superconducting 
neutron star cors". [hep-ph/9910254] 

D.M. Sedrakian, D. Blaschke, K.M. Shahabasyan, D.N. Voskresensky, "Meissner effect 
for color superconducting quark matter". [ hep-ph/0012383 | 

B. Carter, "Application of continuous variation formulas and discrete invariance princi- 
ples to black holes and neutron star models", Annals of Physics, 95, 53 -73 (1974). 

G.A. Vardian, D.M. Sedrakian, "Magnetohydrodynamics of superfluid solutions", Sov. 
Phys. J.E.T.P. 54(5), 919-921 (1981). 



37 



[19] M.A. Alpar, S.A. Langer, J. A. Sauls, "Rapid postglitch spin - up of the superfluid core 
in pulsars", Astroph. J. 282, 533-541 (1984). 

[20] S.E. Thorsett, D.Chakrabarty, "Neutron star mass measurements: radio pulsars", As- 
troph. J. 512, 288-299 (1999). 

[21] E. Chubarian, H. Grigorian, G. Poghosyan, D. Blaschke, "Deconfinement transition in 
rotating compact stars", Astron. Astrophys. 357, 968-976 (2000). lastro-ph/9903489| 

[22] M. Ruderman, "Neutron starquakes and pulsar periods", Nature 223, 597-598 (1969). 

[23] R. Smoluchowski, "Frequency of pulsar starquakes", Phys. Rev. Lett. 24, 923-925 (1970). 

[24] P.W. Anderson, N. Itoh, "Pulsar glitches and restlessness as a hard superfluidity phe- 
nomenon", Nature 256, 25-26 (1975). 

[25] A. Sedrakian, J.M. Cordes, "Spin evolution of pulsars with weakly coupled superfluid 
interiors", Astroph. J. 502, 378-381 (1989). |astro-ph/9802T02 l 

[26] M.A. Alpar, P.W. Anderson, D. Pines, J. Shaham, "Giant glitches and pinned vorticity 
in the Vela and other pulsars", Astroph. J. Lett. 249, L29 (1981). 

[27] A. Sedrakian, I. Wasserman, J.M. Cordes, "Precession of isolated neutron stars. 1. Effects 
of imperfect pinning", Astroph. J. 524, 341-360 (1999). 

[28] M.Q. Alpar, J. A. Sauls, "On the dynamical coupling between the superfluid interior and 
the crust of a neutron star" Astroph. J. 327, 723-725 (1988). 

[29] L. Bildsten, R.I. Epstein, "Superfluid dissipation timescales in neutron star crusts", As- 
troph. J., 342, 951-957 (1989) 

[30] A.D. Sedrakian, D.M. Sedrakian, "Non - stationary dynamics of rotating superfluid sys- 
tems", Sou. Phys. JETP 75, 395-399 (1992). 

[31] M.A. Alpar, P.W. Anderson, D. Pines and J. Shaham, "Vortex creep and the internal 
temperature of neutron stars", Astroph. J. 276, 325-334 (1984); Astroph. J. 278, 791-805 
(1984). 

[32] P.B. Jones, "Rotation of the neutron drip superfluid in pulsars: period discontinuities 
and internal temperatures", Mon. Not. R. Astr. Soc. 246, 315-323 (1990). 

[33] A.D. Sedrakian and D.M. Sedrakian, "Superfluid core rotation in pulsars", Astroph. J. 
447, 305-323 (1995). 

[34] P.B. Jones, "Rotation of the neutron drip superfluid in pulsars: the interaction and 
pinning of vortices" Astroph. J. 372, 208-212 (1991). 

[35] B.K. Link, R.I. Epstein, "Mechanics and energetics of vortex pinning in neutron stars", 
Astroph. J. 373, 592-603 (1991). 

[36] B.K. Link, R.I. Epstein, "Thermally driven neutron star glitches", Astroph. J. 457, 844- 
854 (1996). 



38 



[37] A. Sedrakian, J.M. Cordes, "Vortex - interface interactions and generation of glitches", 
Mon. Not. R. Astr. Soc. 307, 365 (1999). [|astro-ph/980604^ 



[38] M. Ruderman, "Crust breaking by neutron superfluids and the Vela pulsar glitches", 
Astroph. J. 203, 213-222 (1976). 

[39] M. Ruderman, "Neutron star crustal plate tectonics. I. Magnetic dipole evolution in 
millisecond pulsars and low mass X-ray binaries", Astroph. J. 366, 261-269 (1991). 

[40] B. Carter, D. Langlois, D.M. Sedrakian, "Centrifugal buoyancy as a mechanism for neu- 
tron star glitches", Astron. Astroph 361, 795-802 (2000). |astro-ph/0004121| . 



[41] R. Prix, "Slowly rotating two-fluid neutron star models", Astron. Astrophys. 352, 623- 
631 (1999). Hastro-ph/9910293H . 

[42] S.L. Shapiro, S.A. Teukolsky, Black holes, white dwarfs, and neutron stars (Wiley, New 
York, 1883). 

[43] F. Pacini, "Energy emission from a neutron star", Nature, 216, 567-568 (1967). 

[44] P. Goldreich, W.H. Julian, "Pulsar electrodynamics", Astroph. J. 157, 869 (1969). 

[45] N. K. Glendenning, S. Pei, F. Weber, "Signal of quark deconfinement in the timing 
structure of pulsar spin down", Phys. Rev. Lett. 79 , 1603-1606 (1997). 

[46] V.S. Beskin, A.V. Gurevich, Ya. N. Istomin, Physics of the pulsar magneto sphere (C.U.P., 
Cambridge, 1993). 

[47] A. Melatos, "Spin down of an oblique rotator with a current starved outer magneto- 
sphere", Mon. Not. R. Astr. Soc. 288, 1049-1059 (1997). 

[48] A.G. Lyne, R.S. Pritchard, F.G. Smith, "Crab pulsar timing 1982 - 87", Mon. Not. R. 
Astr. Soc. 233, 667-676 (1988). 

[49] A.G. Lyne, R.S. Pritchard, F.G. Smith, F. Camilio, "Very low braking index for the Vela 
pulsar", Nature 381, 497-498 (1996). 

[50] J. Shaham, "Free precession of neutron stars: the role of possible vortex pining" , Astroph. 
J. 214, 251-260 (1977). 



[51] A. Melatos, "Radiative precession of an isolated neuton star". | astro-ph / 000403TI 



[52] B. Carter, H. Quintana, "Relativistic formulation of the neutron starquake theory pulsat 
glitches", Annals of physica 95, 74-89 (1975). 

[53] B. Carter, H. Quintana, "Stationary elastic deformation of a relativistic neutron star 
model", Astroph. J. 202, 511-522 (1975). 

[54] B. Carter, D. Langlois, "Kalb-Ramond coupled vortex fibration model for relativistic 
superfluid dynamics", Nuclear Physics B454, 402-424 (1995). |hep-th/9611082 l 



39 



B. Carter, "Relativistic dynamics of vortex defects in superfluids" in Topological defects 
and the non- equilibrium dynamics of symmetry breaking phase transitions (Les Houches 
99), ed. Y.M. Bunkov and H. Godfrin, pp 267-301 (Kluwer, 2000). fer-gc/990703g l 



D. Langlois, D. M. Sedrakian, B. Carter, "Differential rotation of relativistic superfluid 
in neutron stars" Mon. Not. Roy. Astr. Soc. 297, 1198-1201 (1998). [|astropri/9711042 | 



L. Lindblom, G. Mendel, "The oscillations of superfluid neutron stars", Astroph. J. 421, 
689-704 (1994). 

U. Lee, "Nonradial oscillations of neutron stars with superfluid cores", Astron. Astrop. 
303, 515-525 [1995). 

A. Sedrakian, I. Wasserman, "Perturbations of self - gravitating, ellipsoidal superfluid - 



normal fluid mixture", Phys. Rev. D63, 024016 (2000). lastro-ph/000433l| 



G.L. Comer, D. Langlois, L.M. Lin, "Quasi-normal modes of general relativistic neutron 
stars", Phys. Rev. D60, 104025 (1999). fer-qc/ 9908040 ] 



N. Andersson, G.L. Comer, "Slowly rotating general relativistic superfluid neutron stars" . 
fer-qc/00090"89[ 



N. Stergioulas, J. Friedman, "Non axisymmetric neutral modes in rotating relativistic 
stars", Astroph. J. 492, 301-322 (1998). [|gr-qc/9705056 l 



S. Bonazzola, J. Frieben, E. Gourgoulhon, "Spontaneous symmetry breaking of rapidly 
rotating stars in general relativity", Astron. Astrophys. 331, 280-290 (1998). [gr 



qc/9710121 



S. Morsink, N. Stergioulas, S.R. Blattnig, "Quasi - normal modes of rotating relativistic 
stars - neutral modes for realistic equations of state", Astroph. J. 510, 854-861 (1999). 
fer-qc/98060"08 | 



E. Gourgoulhon, P. Haensel, R. Levine, E. Paluch, S. Bonazzola, J. A. Marck, "Fast 
rotation of strange stars", Astron. Astroph. 349, 851-862 (1999). lastro-ph/9907225 ] 



N. Andersson, D.I. Jones, K.D. Kokkotas, N. Stergioulas, "R-mode runaway and rapidly 
rotating neutron stars", Astroph. J. 534, L75 (2000). [|astro-ph/0002114 l 



G.E. Volovik, "Exotic properties of superfluid helium 3" (World Scientific, Singapore, 
1992). 

J. A. Sauls "Superfluidity in the interior of neutron stars" in Timing Neutron Stars, Nato 
ASI C262, ed. H. Ogelman, E.J. P. van den Heuvel. pp 457-490 (Kluwer, Dordrecht, 
1988). 

G. Mendel, L. Lindblom, "The coupling of charged superfluid mixtures to the electromg- 
netic field", Annals of Physics 205, 110-129 (1991). 

B. Carter, D. Langlois, "Relativistic models for superconducting superfluid mixtures", 
Nuclear Phys. B531, 478-504 (1998). [|gr-qc/980602l 



40 



[71] B. Carter, "Superconducting superfluids in neutron stars", Gravitation and Cosmology 
6, Supplement, 204-213 (2000). [|astro-ph/OOloT09l 



[72] I.L. Bekarevich, I.M. Khalatnikov, "Phenomenological derivation of vortex motion in 
Hell", Sov. Phys. J.E.T.P. 13, 643-646 (1961). 

[73] B. Carter, "Covariant Theory of Conductivity in Ideal Fluid or Solid Media" , in Relativis- 
ts Fluid Dynamics (C.I.M.E., Noto, May 1987) ed. A.M. Anile, & Y. Choquet-Bruhat, 
Lecture Notes in Mathematics 1385 pp 1-64 (Springer - Verlag, Heidelberg, 1989). 

[74] V. V. Lebedev, I.M. Khalatnikov, "Relativistic hydrodynamics of superfluids" , Sov. Phys. 
J.E.T.P. 56, 923-930 (1982). 

[75] B. Carter, I.M. Khalatnikov, "Equivalence of convective and potential derivations of 
covariant superfluid dynamics" , Phys. Rev. D45, 4536-4544 (1992). 

[76] B. Carter, D. Langlois, "The Equation of state for cool relativistic two-constituent su- 
perfluid dynamics", Phys. Rev. D51, 5855-5864 (1995). [hep-th/9507 058] 

[77] B. Carter, "Perfect fluid and magnetic field conservation laws in the theory of black hole 
accretion rings", B. Carter in Active Galactic Nuclei, ed. C. Hazard & S. Mitton, pp 
273-300 (Cambridge U.P., 1979). 

[78] B. Carter, "Black hole equilibrium states", in Black Holes (proc. 1972 Les Houches Sum- 
mer School), ed. B. & C. DeWitt, pp 57-210 (Gordon and Breach, New York, 1973). 

[79] W. Unruh, "Experimental black hole evapouration" , Phys. Rev. Letters 46, 1351-1357 
(1981). 

[80] W. Unruh, "Dumb holes and the effects of high frequencies on black hole evapouration", 
Phys. Rev. D51, 2827-2838 (1995). [fer-qc/9409008l 

[81] C. Barcelo, S. Liberati, M. Visser, "Analog gravity from Bose-Einstein condensates". 
fer-qc/0011026 | 



[82] A. Lichnerowicz. Relativistic Hydrodynamics and Magnetohydrodynamics (Benjamin, 
New York, 1967). 

[83] A.H. Taub, "General relativistic variation principle for perfect fluids", Phys. Rev. 94, 
1469-1470 (1954). 

[84] B. Schutz, "Perfect fluids in general relativty: velocity potentials and a variational prin- 
ciple", Phys. Rev. D2, 2762-2773 (1970). 

[85] B. Carter, "Axionic Vorticity Variational formulation for Relativistic Perfect Fluids", 
Class. Quantum Grav. 11, 2013-2030 (1994). 

[86] B.Carter, in A Random Walk in Relartivity and Cosmology, Proc. Vadya - Raychaudhuri 
Festschrift, IAGRG 1983, ed. N. Dadhich, J. Krishna Rao, J.V. Narlikar, C.V. Vishvesh- 
wara, pp 48 -62 (Wiley Eastern, Bombay, 1985). 



41 



[87] B. Carter, I.M. Khalatnikov, "Momentum, Vorticity, and Helicity in Covariant Superfluid 
Dynamics", Ann. Phys. 219, 243-265 (1992). 

[88] W. Israel, "Covariant superfluid dnamics", Physics Letters A86, 79-81 (1981); "Equiv- 
alence of two theories of relativistic superfluid dynamics", Physics Letters A92, 77-78 
(1982) 

[89] W.G. Dixon, Arch. Rat. Mech. Anal. 80, 159 (1982). 

[90] T.S. Olsen, "Instabilities in the Israel-Dixon relativistic superfluid theory", Physics Let- 
ters A149, 71-75 (1990). 

[91] L.D. Landau, E.M. Lifshitz, Statistical Physics (trad. E. and R.F. Peierls), Section VI 
(Pergamon, Oxford, 1959). 

[92] R.I. Epstein, "Acoustic properties of neutron stars" Astroph. J. 333, 880-894 (1988). 

[93] S. Tsuruta, "Thermal properties and detectability of neutron stars", Phys. Rep. 56, 
237-278 (1979). 

[94] A.F. Andreev, E.P. Bashkin, "Three veloity hydrodynamics of superfluid solutions", Sov. 
Phys. J.E.T.P. 42, 164-167 (1976). 

[95] O. Sjoberg, "On the Landau effective mass in asymmetric nuclear matter", Nucl. Phys. 
A265, 511-516 (1976). 

[96] P. B. Jones, "Rotation of the neutron drip superfluid in pulsars: the Kelvin phonon 
contribution to dissipation", Mon. Not. R. Astr. Soc, 257, 501-506 (1992). 

[97] B. Link, R.I. Epstein, G. Baym, "Superfluid vortex creep and rotational dynamics of 
neutron stars", Astroph. J., 403, 285-302, (1993). 

[98] P. Haensel, "Non-equilibrium neutrino emissivities and opacities of neutron star matter" , 
Astron. Astroph., 262, 131-137 (1992). 



42 



